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Abstract 
We  are  dealing  with  a  system  of  two  partial  differential 
equations  of  the  form 

u,  +  f   -  0, 
(1)  ^    "" 

^t  +  Sx  =  0- 
f  and  g  are  functions  of  u  and  v   subject  to  the  following 
conditions : 

(a)  The  system  (l)  is  hyperbolic. 

(b)  The  system  (l)  is  genuinely  nonlinear  in  the  sense  of  [4] 

(c)  The  interaction  of  two  shocks  of  the  same  family  always 
produces  a  shock  of  the  same  family  plus  a  rarefaction  wave  of  the 
opposite  family. 

Denote  by   e   the  oscillation  of  the  initial  data.   The 
initial  data  need  not  have  bounded  total  variation.   Our  results 
are  as  follows: 

Theorem  5 • !•   For  epsilon  small  enough  the  solution  of  the 
initial  value  problem  exists  for  all  positive  time. 

Theorem  5 -2 •   If  the  initial  data  are  periodic  the  solution 
decays  like   l/t . 

The  proofs  are  based  on  methods  introduced  in  [2]  supple- 
mented by  further  inequalities  which  take  into  account  the  inter- 
action of  shock  waves  and  rarefaction  waves  of  the  same  family. 
The  new  estimates  prove  that  such  interactions  must  occur  and 
that  when  they  do  occur  they  lead  to  decay  in  the  solution. 
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The  inequalities  presented  in  this  paper  can  also  be 

used  to  prove  if  the  initial  data  are  constant  outside  of  a 

1/2 
compact  interval,  the  solution  decays  like   l/(t) 
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DECAY  OF  SOLUTIONS  OF  SYSTEMS  OF 
HYPERBOLIC  CONSERVATION  LAWS . ^ 


by 

James  Glimm'"  and  Peter  D.  Lax 


Preface 
In  this  paper  we  study  coupled  systems  of  nonlinear 
conservation  equations: 

u^  4  f^  =  0,   v^  +  g^  =  0, 

f,    g   functions  of  u  and  v.   Me  assume  that  this  system 
is  strictly  hyperbolic,  i.e.  has  two  real  and  distinct 
characteristic  speeds   p   and  X'      These  characteristic 
speeds  depend  on  u   and  v;   it  is  well  known  that  this 
causes  in  general  the  appearance  of  discontinuities,  i.e. 
shock  waves.   Our  aim  is  to  show  that  the  dependence  of  the 
characteristic  speeds  on  u  and  v  also  causes  shock  waves 
and  rarefaction  waves  of  the  same  family  to  r\an  into  each 
other  and  to  cancel  each  other  out.   Our  main  estimate. 
Theorem  4.1,  is  a  quantitative  form  of  this  observation. 
This  estimate  is  used  to  prove  Theorem  ^.1,    the  existence  of 
solutions  with  arbitrarily  prescribed  bounded,  measurable 
initial-data,  with  small  oscillation.   Another  application 

The  results  of  this  paper  were  announced  in  [ 3  ] • 
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Is  Theorem  5-2,  accordJng  to  which  space-periodic  solutions 
satisfy  the  decay  law 

I   —  1  _,  const    I   —  I  ^  const 
|u-ul  ^  — :g —  ,  lv-v|  ^  — ^ — 

where   u,   v   denote  the  mean  value  of  u   and  v  per  period. 

Thus  we  find  for  these  equations  that  a  decay  phenomenon 
for  the  solutions  leads  to  a  new  estimate  for  the  solution, 
an  improved  existence  theorem,  and  a  bound  on  the  asymptotic 
behavior  of  the  solutions   for  large  times. 

An  important  ingredient  of  our  method  is  an  analysis  of 
the  interaction  of  the  two  characteristic  modes  of  propagation; 
we  show  that  if  the  solution  differs  little  from  a  constant, 
the  interaction  is  negligible.   Since  this  result  is  used 
essentially  in  our  estimates,  we  can  only  deal  with  weak  signals, 
i.e.  solutions  whose  initial  state  differ  little  from  a  constant. 

¥e  give  now  an  intuitive  description  of  the  ideas  involved 
illustrating  them  first  on  a  single  conservation  equation 

(P.l)  ^t  "^  ^x  "  °' 

The  differentiated  form  of  the  equation  is 

u,  +  f '(u)u^  =  0. 

We  take  for  f  the  simplest  nonlinear  function: 

f(u)  =  p-  u  ,    a  positive  constant. 
The  nonconservation  form  of  (P.l)  is  then 
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( P.2 )  u,  +  auu   =0, 

which  shows  that  u   is  constant  along  curves  satisfying 

(P. 3)  ^  -  au. 

Since   u   is  constant  along  them,  these  curves,  called 
characteristics,  are  straight  lines.   In  general  these  straight 
lines  intersect;  at  a  point  of  intersection  the  solution  would 
be  multivalued. 

Instead  of  becoming  multivalued  the  solution  becomes 
discontinuous,  satisfying  equation  (P.l)  in  the  sense  of 
distributions.   For  piecewise  continuous  solutions  this  means 
the  fulfillment  of  the  jump  relation 

(P.4)  ^=IlI  =  s 

^*   [u] 

where  x  =  x{t)   is  the  curve  of  discontinuity,  and 
[f]  =  fp  -  f^  >       [u]  =  u^  -  u^ 

is  the  difference  between  the  values  of   f,   respectively  u, 
on  the  right  and  left  side  of  x(t).   For  our  special  choice 
of   f   (P. 4)   becomes 

(P  5)  ^  -  a  ""-^^""^ 

It  is  well  known  that  only  those  discontinuous  solutions 
are  acceptable,  either  mathematically  or  physically,  which 
satisfy  the  so-called  entropy  condition,  which  requires  that 


IX 


(P. 6) 


f'(uj  >  lil  >  f'(u^) 


'i 


[u] 


Geometrically  this  condition  asserts  that  the  characteristics 
on  either  side  of  the  discontinuity,  when  drawn  in  the  direction 
of  increasing   t,   run  into  the  curve  of  discontinuity: 


t/\ 


Shock. 


p-  au  ,  the  entropy 


Fig.  P.l 

Discontinuities  satisfying  (P. 6)  are  called  shocks. 

Note  that  for  our  choice  of   f 
condition  takes  the  form 


u  .+u 
au^  >  a  — 2 —  ^  ^^r' 


which  for  a  >  0   is  equivalent  with 


(P.?) 


Ufl  >  u  . 
Ji  r 


We  return  now  to  continuous  solutions;  let  x-,(t)   and 
Xp(t)   be  a  pair  of  characteristics,   0  ^  t  <_  T;   it  follows 
from  the  description  of  u  as  constant  along  characteristics 
that  u(x,0)   on  the  interval   [x-,  ( O)  ,Xp(  O)  ]   and  u(x,T)   on 
the  interval   [x-,(T),  Xp(T)]   are  equivariant .   In  particular 


X 


the  increasing  and  decreasing  variation  respectively  of  u   on 
these  two  Intervals  are  the  same. 

More  generally,  let  <r       and   oi   be  any  pair  of  curves 
connecting  x-,(t)   and  Xp{t),   neither  of  them  tangent  to  a 
characteristic;  then  likewise  u   on   or:   is  equivariant  with 
Op  . 

Suppose  that  there  are  shocks  present;  let   y-i(t)   and 


u   on   o^ 


y^ ( t )   be  a  pair  of  characteristics  which  run  into  shocks  at 


times   t 


-,   and   tp,   respectively.   Define   y-(t)   for  t  >  t. 


as  the  shock  into  which  y.   has  run.   Let  T  be  some  value 

1 

>  t.;   there  are  two  points   x-,   and  Xp   in  the  Interval 
[y-i  (  O)  jYoI  O)  ]   such  that  the  characteristics  issuing  from  x. 
intersects  the  shock  y.   at  time   T: 


T   .. 

t^ 


y 


x^(T)  =  y^(T) 


1 


X- 


Xo  y, 


X 


Fig.  P. 2 
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If  there  are  no  further  shocks  contained  between  x-,   and  Xp 
for  0  ^  t  _<  T   then  u(x,0)   on  the  open  intervals   (x-,(0), 
Xp(0))   and  u(x,T)   on  the  open  interval  (x-,  (T)  ,Xp(  T)  )   are 
equivariant  ;  otherwise  u(x,T)   has  a  smaller  variation. 
This  shows  that  the  effect  of  characteristics  running  into 
shocks  is  to  decrease  the  total  variation  u(x,t).   ¥e  make 
now  a  quantitative  estimate  of  this  decrease. 

Denote  by  D(t)   the  width  of  the  strip  bounded  by 
X-,  ( t )   and  Xp(  t )  : 

D(t)  =  x^{t)    -  x^(t). 

Differentiating  with  respect  to  t  we  get,   upon  using  relation 
(P.3),   that 

dt  °^*^  "  dt  ^2  "  dt  ^1  "  ^^2  -  ^^1- 

Since  Up  and  u-,   are  constant  along  characteristics,  we 
deduce  after  integrating  with  respect  to  t   that 

(P. 8)       D(T)  =  D(0)  +  a(u2-u^)T. 

Let  L   be  any  interval  on  the  x-axis;  we  subdivide  it 
into  subintervals  [y  ■_-,  fY  ■]  ,    j=l,  .  .  .n,  on  which  u(x,0)   is 
monotonic,  increasing  for   j   odd,  decreasing  for   j   even. 
As  before  we  define   y.(t)   as  the  characteristic  issuing  from 
y.,   continued  as  a  shock  if  it  runs  into  a  shock.   Denote  by 
L(t)   the  interval   [y^l t ) ,y^( t ) ] . 

The  analysis  presented  before  shows  that  for  all   t 
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u(x,t)   is  increasing  on   (  y  •_-,  ( t )  ,y  .(t )  )   for   j   odd,  decreasing 

for   j   even;  according  to  (P.7)   u   decreases  across  shocks. 

Therefore  the  total  increasing  variation  of  u(x,t)   across 

L(t)   is 

(P. 9)  H   u  (t)  -  u.  .(t), 

j  odd   ^       ^ 

where  ^i-_-|(t)  denotes  the  value  of  u   on  the  right  edge  of 

y._-|(t),   u.(t)  the  value  of  u   on  the  left  edge  of  y.(t). 

Denote  as  before  by  ^.-.-[(t),  x.(t)   characteristic 

starting  inside  y  •_-]  ^  y-   which  intersect   y.  -.(t),   respectively 

y.(t)   at   t  =  T.   u.(t)   is  then  the  constant  value  of  u   on 

J  J 

x.(t).   Denoting  x.(t)  -  x._-,(t)   by  D.(t)   we  obtain  from 
(P.8)  that 

D.(T)  =  D.(0)  +  a[u.(t)-u.   (t)]T 

Summing  over  odd   j   we  get 

XI  Dj(T)  =  Y^   Dj(0)  +  aT  ^  u^.  -  u._^. 

The  Intervals   [x.,x.  ,]   are  disjoint   and  lie  in  L(t): 
therefore  the  sum  of  their  length  cannot  exceed  the  length 
|l(T)1.   Using  this  and  (P. 9)  we  can  turn  (P. 10)  into  the 
following  inequality: 

(P. 11)     Total  increasing  variation  of  u(x,T)   over  L(T) 

<  |L(T)|  . 
aT 

Suppose  in  particular  that  u(x,t)   Is  periodic  in  x  with 
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period  p;   take   L   to  be  of  length  p;   then  by  periodicity 
L(t)   has  length  p   for  all   t   and  we  deduce  from  (P. 11)  that 

/p  -in\  Total  variation  of  u(x,t)   per  period    1 

i.e.   that  periodic  solutions  decay  like   1/T.   It  is  remarkable 
that  the  bound  on  the  right  side  is  independent  of  the  particular 
solution  considered. 

¥e  wish  to  generalize  our  method  to  systems  of  two  equations. 
It  is  well  known  that  there  exist  Riemann  invariants   w  and   z 
which  are  functions  of  u  and  v  and  which  satisfy  equations 
of  the  form 

w^  +  pw^  =  0,   ^t  "*"  '^^x  "  ^' 

where   p   and   z   are  functions  of  w  and   z.   These  equations 

express  the  constancy  of  w,   respectively  z,      along  curves 

x(t)   satisfying 

dx  ,  .   -,    dx   -, 

^  =  p  ,   respectively  -Tf-  =  ^' 

These  curves  are  called   p-characteristics  and   A-characteristics . 
It  follows  as  before  that  if  x^(t),  Xp(t)   are  two  p-character- 
istics, then  u(x,0)   along   [x-,  ( O)  ,Xp(  O)  ]   and  u(x,T)   along 
[x-,  (T)  ,Xp(T)  ]   are  equivariant. 

The  p-characteristics  are  no  longer  straight  lines;  never- 
theless it  can  be  shown  that  in  general  they  intersect,  provided 

that   p  genuinely  depends  on  w,  i.e.  that   p  j^  0.      This  leads 

w 

to  the  formation  of  shock  waves.   The  analogue  of  condition  (P. 6) 
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for  systems  is 

(P. 13)     p(u^,v_g)  >  s  >  p(u^,v^), 

where   s   is  the  speed  of  propagation  of  the  discontinuity. 
Discontinuities  satisfying  (P.I3)  are  called   p-shocks; 
similarly  there  are   A-shocks.   The  geometric  interpretation 
of  (P. 13)  is  that   p-characteristics  drawn  in  the  direction  of 
increasing  t   run  into   p-shocks;  from  this  we  deduce  as  before 
that  the  presence  of   p-shocks  decreases  the  total  variation 
of  w. 

To  proceed  further  we  have  to  estimate  also  the  effect 
of  A-shocks  on  the  total  variation  of  w.   Here  it  is  important 
to  recall  that  across  weak  A-shocks  the  change  in  w   is  pro- 
portional  to  (change  in  z)^ ;      this  is  basic  to  the  following 
estimate 

(P.l4)     Change  in  total  variation  of  w   due  to  A-shocks 

^  0(e)V^ 

where   e   is  the  oscillation  of  the  solution,   V   its  initial 

total  variation. 

There  is  a  second  difficulty  we  have  to  cope  with: 

The  differential  equation  for  the  width  of  a  strip  bounded 

by   p   characteristic,   D(t)  =  yi^i"^)    -   x-,(t),   is  as  follows: 

4  D  =  P2  -  P;i_  =  P^(w2-w-^)  +  p^(z2-z-^), 
where  on  the  right  we  have  applied  the  mean  value  theorem.   Since 
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we  are  dealing  with  solutions  whose  total  oscillation   e   is 
small,   p   and   p   themselves  oscillate   0(e).   The  quantities 
w.   and  Wp   are  independent  of   t,   but   z-,   and   Zp   are  not  I 
We  can  overcome  this  difficulty  by  measuring  the  width  of  the 
strip  bounded  by  p-characteristics   not  between  points  with  the 
same   t   coordinates  but  between  points  which  lie  approximately 
on  the  same  A-characteristics .   Since   z   is  constant  along 
A-characteristics  this  makes  the  quantity   Zp  -  z-,   appearing 
in  the  differential  equation  small. 

In  Section  1  we  carry  out  the  details  of  estimating  D(T) 
for  solutions  without  shocks. 

In  Section  2  we  construct  approximate  characteristics 
for  difference  approximations  introduced  in  [ 2  ] .   We  derive 
approximate  conservation  laws  of  the  increasing  and  decreasing 
variation  of  w  and   z;   these  laws  are  formulated  as  a  balance 
between  the  amount  of  shock  wave  and  rarefaction  wave  of  either 
kind  entering  and  leaving  a  region,  the  amount  of  rarefaction 
and  shock  wave  of  the  same  family  cancelling  each  other  in  the 
region,  and  a  correction  term  accounting  for  the  interaction 
between  waves  belonging  to  different  families . 

In  Section  5  we  construct  by  passage  to  the  limit  character- 
istic curves  and  shock  curves  for  exact  solutions,  and  derive 
inequalities  for  the  variations  of  w  and   z. 

In  Section  4  we  study  the  spreading  of  strips  bounded  by 
characteristics  in  the  presence  of  shocks  and  deduce  an  analogue 
of  inequality  (P.ll). 
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In  Section  5  we  use  this  Inequality  to  construct  solutions 
whose  initial  values  are  not  of  hounded  total  variation.   For 
periodic  solutions  we  deduce  the  analogue  of  (P. 12). 

A  word  of  caution:  the  proof  of  lemma  3.4  is  heavy  going 
and  we  advise  the  reader  to  read  the  rest  of  this  paper  first. 

As  we  have  indicated,  we  can  only  deal  with  solutions 
whose  oscillation  is  small.   Recently,  [5],  Nishida  succeeded 
in  solving  the  initial  value  problem  in  the  large  for  the 
special  system:   f  =  v,  g  =  -a/u.   It  would  be  interesting  to 
study  the  decay  of  solutions  of  the  system  considered  by 
Nishida. 
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1.   The  Widening  of  Rarefaction  Waves. 

In  this  section  we  carry  out  rigorously  the  idea  alluded 
to  at  the  end  of  the  preface:  we  obtain  a  lower  ho-und  for  the 
distance  of  two  characteristics   In  terms  of  the  amount  of 
rarefaction  wave  contained  between  them;  here  we  only  treat 
shockless  solutions. 

The  equations  under  consideration  are  of  the  form 
(1.1)  u^  +  f^  =  0,   v^  +  Sx  =  0 

where   f  and  g   are  nonlinear  functions  of  u  and  v. 
(l.l)   can  be  written  as  a  quasillnear  system 

u,  +f  u   +f  V   =0, 

t     U   X     V   X     ' 

(1.2) 

^t  +  Su  ^x  +  S^  ^x  =  °- 

We  assume  that  (1.2)  is  hyperbolic,  i.e.  that  the  coefficient 
matrix  of  the  x-derivatives. 


■'f      f 

U       V 


=u    °v 


has  real  and  distinct  eigenvalues  in  some  neighborhood  of  a 
point  which  we  take  to  be  the  origin   (0,0).   It  is  well  known 
that  for  any  quasilinear  hyperbolic  system( 1 . 2)for  two  functions 
Riemann  invariants  can  be  introduced  as  new  variables;  call  these 
w  and   z.   In  terms  of  these  the  differential  equations  (1.2) 
take  on  the  simple  form 

w,  +  p  w  =  0 

(1.3) 

z^  +  A  z^  =  0 

where   p   and  A   are  functions  of  z   and  w.   We  assume  that 

p   and   A   genuinely  depend  on  w  and   z   respectively,  i.e. 

that   p   and   A   are  nonzero.   Replacing  w  by   -w,  z   by 
"^  w        z 

-z   if  necessary  these  conditions  of  genuine  nonlinearity  can 
be  put  in  the  form 

(1.4)  p^  >  0,   A^  >  0. 

A   and   p   are  called  characteristic  speeds;  curves   x(t) 
propagating  with  characteristic  speed,  i.e.  satisfying 

,-,    r-\  dx   ^      dx 

(1.5)  —  =  A   ,   =  p 

dt         dt 

are  called  A-,   respectively  p-characteristics .   Since   A 
and   p   vary  from  solution  to  solution,  so  do  characteristic 
curves . 

We  assume  that   A  and   p   are  distinct,  say  A  "^  p .  Since 
in  what  follows   w  and   z   differ  little  from   (0,0), 
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we  may  without  loss  of  generality  take   A  <  0  <  p;   then 
all  ^-characteristics  move  to  the  left,   p-characteristics 
to  the  right;  this  contributes  to  the  clarity  of  diagrams. 

Equation  (1.3) states  that   w   is  constant  along   p-charac- 
teristics,  z   constant  along  A-characteristics .   This  fact 
is  basic  to  all  that  follows;  we  shall  use  it  to  prove  the 
following 

Theorem  l.l;   Suppose  that   A   and   p   are  twice  dif ferentiable 
functions  of  w  and   z.   Let   z(x,t),   w(x,t)   be  a  differentiable 
solution  of(1.3)which  is  periodic  in   x  with  period  p   and  which 
exists  for   0  ^  t  j<  T .   Then . 

T  <  ^Q^^^  P  +  p 

~   CD  +  C 

where 

CD  =  total  variation  per  period  of  w(x,0) 

^  =  total  variation  "     "     "   z(x,0) 
provided  that 

e  =  Max^|w(x,0)  I,  |z(x,0)  K 

is  less  than  some  constant. 

Remark:   Since  solutions  are  uniquely  determined  by  their 
initial  data  it  suffices  to  assume  that  the  initial  values  of 
w  and   z   are  periodic. 

Proof;   From  the  constancy  of  Riemann  invariants  along 
characteristics  it  follows  from  the  definition  of  e   that 

(1-6)  |z(x,t) I  <  e,  |w(x,t) I  <  e 

for  all   X  and  all   t. 


Let  X-,   and  Xp   be  any  two  points,  say  x,  "^  x„, 
denote  by  x-,(t)   and  Xp(t)   the   p-characteristics  with 
respect  to  the  solution  under  consideration  issuing  from  these 
points : 


Fig.  1.1 


Denote   A(0,0)   by  A  :   we  draw  a  spacelike  straight  line 
propagating  with  a  speed  A-,   through  the  point   (xp(t),t)   and 
denote  by  t*(t)   the  time  at  which  this  line  intersects  the 
curve  x-,(t),   see  Fig.  1.   t*   satisfies 


(1.7) 


X2(t)  -  x^(t*)  =  A^(t-t*) 


Differentiating  with  respect  to   t  we  get 


(1.8) 


P2  -  Pi-t*  -  A-^d-t-^) 
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where   •   denotes  differentiation  with  respect  to   t   and   p, 

•  « 

and   pp   abbreviate   x-,   and  Xp .   Solving  for   t'^   we  get 

(1.9)  t*  =  — — -  =  1  +  — — i  . 

Pl-^1       Pl-^1 

We  denote  by  D*(t)   the  horizontal  width  of  the  strip  between 
x^  ( t* )   and  Xp{t)   in  the   A-,   direction,  i.e. 

(1.10)  D*(t)  =  X2(t)  -x-^(t*). 

Differentiating  (I.IO)  with  respect  to   t  we  get,  using  (1.9)j 
that 

(1.11)  ^  ^*  =  P2  -  Pi*""  =  (P2-Pl^  —  =   tP^  — 

\=Pj_  \-pi 

Here   [p]   abbreviates   pol^)  -  p-|(t*). 
We  choose  now 

A-,  =  A(0,0)  -  const  e, 

the  constant  chosen  so  large  that  the  line  with  slope   A-,   is 
spacelike.   This  is  possible  since   |w|   and   |z|   never  exceed 
e;      for  the  same  reason,   p.,  =  p(0,0)  +  0(e).   Using  the 
abbreviations 

p(0,0)  =  p^,      A(0,0)  =  Aq 


we  have  therefore 

A.       A 


o 


^I'Pl    ^o-P( 
so  (1.11)  can  be  rewritten  es 


(l+0(e)); 


(1.12)  ^  ^*  =  tp]  x-^  (l+0(e)) 

o  Po 


(1-15) 
we  get 


Using  the  mean  value  theorem  and  the  abbreviations 
[w]  =  w(x2(t),t)  -  w(x-^(t*),t*) 
[z]  =  z(x^(t),t)  -  z(x.(t*),t*) 


[p]  =  a[w]  +  b[z] 


where  a  =  p  ,   b  =  p   at  some  intermediate  points.   We 
substitute  this  into  (1.12);   we  can  write  the  result  as 

(1.14)  ^  D*  =  c^[w]{l+0(e))  +  g[z] 

where 

Cq  =  a^  ^  ,  g  =  b   °   =  b  2_  (i+0(e)). 

^o-Po         ^1-Pl     ^Q-Pl 

Since  according  to  (1.4)   p,^   is  positive,  so  is   c^. 

Next  we  estimate  [z];   for  this  we  make  use  of  the  constancy 
of  z  along  A-characteristics : 

(1.15)  z(x^(t^),t*)  =  z(x2(t-tQ),t-tQ) 

where   t  -  t    is  the  time  when  the  backward   A-characteristic 
through  x-,(t*),t*   intersects  the  curve  Xp(t),   as  shown  on 
figure  1.2  below: 


X2(t-t^) 


Fig.  1.2 


Since  w  and  z   oscillate  at  most   e,   the  speed  A 
differs  at  most  by  0(e)   from  A,   chosen  earlier.   The  time 
difference  t   is  hounded  by  the  difference  of  the  reciprocals 
of  the  speed  times  the  distance  to  be  travelled: 


(1.16) 


tQ(t)  =  0(e)D*(t) 


using  (1.13)  and  (I.I5)  we  can  write 


[z]  =  z(x2(t),t)  -  z(x-j_(t*),t*) 


=  z(x2(t),t)  -  z(x2( t-t^,t-t  ) 

"  /    ds  z(x2(s),s)ds 
t-t 


Substituting  this  into  (l.l4)  gives 
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dt 


D*(t)  =  c[w]  +  &     I      z^(s)ds 


t-t 


where 


and 


^  =  d?  ^(^2^^)'^) 


c  =  c^(l+0(e)) 


Since  w   is  constant  along  both  characteristics   x-,(t)   and 

Xp(t)   it  follows  that  so  is   [w] .   Therefore  integrating  the 

above  with  respect  to   t   from  0  to  T  gives 

T  t 
(1.17)      D*(T)  =  D*(0)  +  c[w]T  +  /  /  |gZg(s)|dsdt. 


0  t-t 


o 


To  estimate  the  last  term  on  the  right,  we  interchange 

the  order  of  integration 

T  t  T  S+t;j_(  s) 

(1.18)      ff     |g(Zg)|(s)dsdt  -jTJ      |g|  |z^|dt  ds  +  q 

6  t-t  0  's 

o 

where   t-,   is  defined  in  figure  3  below: 


and   q   denotes  the  Integral  over  a  triangular  piece: 

0    s+t2_(s) 
(1-19)  q  =f        f  Igz^ldt  ds. 

-%{0)    0 

To  estimate  the  right  side  of  (l.l8)  we  need  an  estimate  for 
t-,;   from  figure  1-3  we  get  the  relation 

0  ^  t-^(s)  =  t^(s+t-^(s)). 

Using  (l.l6)  we  get 

t^(s)  =  0(e)D*(s+t-^(s)). 

Since  D*   is  Lipschitz  continuous 

D*(s+t^(s))  =D*(s)  +  0(t^(s)); 
this  combined  with  the  above  estimate  for  t-,   gives 

t^(s)  =  0(e)D*(s). 

The  function  g   is  uniformly  bounded;  therefore  using  the  above 

estimate  for  t-,   gives  the  following  upper  bound  for  the  double 

integral  on  the  right  side  of  (l.l8): 

T 
0(e)j  |Zg(s)  |D*(s)ds. 
6 

We  can  estimate   q   in   (1.19)  similarly: 

(1.20)  |q|  ^  0(e)D*(0)V 

where 


0 


(1.21)  V  -  y   Izgl'^s 

-t;(o) 

Combining  (1.17)  and  (l.l8)  and  using  the  above  estimates 

we  get  the  following  two-sided  estimate  for  D*(T): 

T 
[l-O(e)V)] [D*(0)  +  c[w]T]  -  0( e ) j  1 Zg( s ) |d*( s ) | ds  <  D*( T) 

(1.22)  °    ^ 

^  [l+0(e)V)] [D*(0)  +  c[w]T]  +  0(  e ) T  | Zg( s ) |d* ( s)ds . 

6 

Next  we  make  use  of  the  following  classical  inequality: 

Lemma:   Let  m  be  an  increasing  function  and  suppose  that 

for  each  positive   T,  M   satisfies  the  inequality 

T 

(1.23)  M(T)  ^  F(T)  +  rM(t)dm(t) 

6 

Then  M   satisfies  the  inequality 

(1.24)  M(T)  <  M-^(T) 

where  M-,   is  the  solution  of  the  integral  equation  obtained 

from  (1.25)  by  replacing  the  sign  of  inequality  by  equality: 

T 
(1.23)-L  M-l(T)  =  F(T)  +  I  M^(t)dm(t)  . 

0 

The  explicit  expression  for  M,   is 

T 

(1.25)  M-^(T)  =  e"^^T)p(o)  ^  J^m(T)-m(t)  dF  ^^. 

6 

here  m(0)   is  taken  to  be   0. 

Suppose  that  in  addition  to  (1.23)   M   satisfies  an  inequality 


10 


in  the  opposition  direction: 

T 
(1.26)       G(T)  -  r  M(t)din  <  M{T)  . 

6 

Substituting  the  upper  "bound  (1.24)  into  the  left  side  of 

(1.26)  and  using  (1.23)-j^  we  obtain  the  following  lower  bound 
for  M: 

(1.27)  G(T)  +  F(T)  -  M^(T)  ^M(T). 

Suppose   F(t)   is  an  increasing  function;  then  the  integral  on 
the  right  in  (I.25)  will  increase  if  the  exponent  m(T)  -  m(t) 
is  replaced  by  the  larger  quantity  m(T),   and  we  get  the 
following  upper  bound  for  M-,  : 

M-l(T)  <  e"^^^^F(T) 

Substituting  this  into  (1.27)  gives 

(1.28)  G(T)  +  (l-e'^^^^)F(T)  <  M(T). 

¥e  apply  this  to  the  function  D*   as  M;   inequality 

(1.22)  is  inequalities  (1.23)  and  (I.26),  with 

t 

(1.29)  m(t)  =  0(e)J  |dz(s)  I  =0(e)Z(t) 

0 
F(T)  =  [l+0(e)V][D*(0)  +  c[w]T] 

G(T)  =  [l-0(e)V] [D*(0)  +  c[w]T]. 
Inequality  (1.28)  in  this  case  gives 
(1.50)       [2-e°^^^^^^^-  0(e)V] [D*(0)+c[w]T]  <  D*(T) 
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where   c  =  [l-0(e)]c  . 

The  quantity   Z(t)   defined  in  (1.29)  is  the  total  variation 
of   z   along  the  characteristic   Xp(t).   Using  again  the  constancy 
of  z  along  A-characteristics  we  deduce  the  following  law  of 
conservation  of  total  variation: 

Let   C-,   and   Cp   be  two  curves,  both  transversal  to 
A-characteristics,  whose  left  endpoints  and  right  endpoints 
respectively  lie  on  the  same  A-charact eristic  (see  figure  1.4): 


A -characteristics 


Then  z   along  C 


1 


and 


Fig.  1.4 


Cp  are  equivariant. 


It  follows  from  this  that  the  total  variation  of   z   along 
C-.   equals  its  total  variation  along  Cp. 

¥e  apply  this  result  to  the  case  when  Cp   is  the 
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p-characteristic   Xp   between   t  =  0   and   t  =  T,   and   C-, 
is  the  Interval  of  the  x-axis  "between  Xp   and   y{T)   defined 
as  that  point  where  the   A-characterlstlc  through  the  point 
Xp(T),T   Intersects  the  x  axis,  see  figure  1-5: 


X. 


y(T) 


Fig.  1.5 

Clearly  the  length  of  the  Interval   [Xp,y(T)]   Is   0(T). 

We  recall  now  that  the  total  variation  of  the  initial 
value  of  z  per  period  is  denoted  by  C;   therefore  the  total 
variation  of  z   along  an  interval  of  length  L  is  L(!^/p. 
Thus  for  T  >  0(p) 

Z(T)  =  Total  variation  of  z   along  Xp(t)  = 
total  variation  of   z   in   (Xp,y)  =  0(CT)/p• 
A  similar  estimate  for  the  quantity  V  defined  in  (1.21) 
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gives  V  _^  (^   provided  x-,   and  Xp   lie  within  a  single  period. 
Using  these  estimates  for  Z   and  V   in  inequality  (1-30)  gives 

(lOD       [2-e°^^^^^/P-0(e)C]c[w]T  ^  D*(T). 

We  restrict  e,    C   and  T   so  that  this  inequality  is  non- 
trivial,  i.e.  that  the  coefficient  on  the  left  is  positive; 
this  will  be  the  case  if  the  exponent   0(eCT)/p   is  less  than, 
say,  log  1.^,      and   0(e)C   is  less  than  .2: 

(1.32)  T  <  £2Il£i  p   ,   eC  <  const. 

For  T,  e,  C   satisfying  (1.32)  inequality  (I.3I)  implies 

(1.33)  c[w]T  ^  2D*(T) . 

So  far  X-,   and  Xp  were  any  two  points  within  one  period. 
¥e  choose  now  a  set  of  intervals   I"^  =  {x\^    ,xX^^    )    all  within 
one  period,  say   (0,p),  so  that 

i)   the  quantities   [w]'^*  =  w(x^^'  ,0)  -  w'(x^"^'  ,0) 
are  all  positive 

ii)  co/3  1  Yl   [w^*^' 

where,  as  before,   cd  denotes  the  total  variation  of  w  per 
period. 

Such  a  choice  of  intervals  is  possible  since  the  variation 
of  any  function  in  the  positive  direction  is  exactly  one  half 
its  total  variation. 
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The  points   ^-i  >  ^n        were  chosen  to  lie  in  an  interval 
of  length  p.   Consider  the   p-characteristic  curves   x-,'-^(t) 
and  Xp  (t)   issuing  from  these  points.   Since  the  solution 
under  consideration  is  periodic  with  period   p,   it  follows 
that  for  any  t   the  points   x-.'^lt)   and   Xp  (t)   also  lie  in 
an  interval  of  length   p.   Denote  by  D.(t)   the  width  of  the 

J 

strip  bounded  by  the   j    pair  of  characteristics: 

Dj.(t)  =  Xg-^'lt)  -  x-^^-^^t). 
It  follows  from  the  above  analysis  that 

^-D.{t)    <   p. 

We  recall  now  the  quantity  D.*{t)   which  measures  the  width 

J 

of  the   j    strip  in  the  A-,   direction.   Clearly 

D  *(t)  ^  const  D.(t) 

so  STomming  we  get   from  the  above   inequality  that 

XI  D.-*(t)   ^  const  p. 
J 

On  the  other  hand  if  T  satisfies  (1.32)  then  inequality  (1.33) 
holds  for  each  strip: 

c[w]'^'t  ^  2D  .*(T)  . 
(J 

Stunming  with  respect  to  j   and  using  the  facts  that 

Y~   D.*(t)  <_  const  p  and  that  ^  [w]  "^  >  03/3  we  get  that 

(1.3^)  a)T  <  const  p 
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provided  that  (1-32)  is  satisfied.   Dividing  by  cd  we  get 
(1.35)  T^const_p^ 

Since  the  roles  of  w  and   z   are  interchangeable,  we  may 
assume  that   |  _<  cd.   In  this  case  for   e   less  than  a  constant 
independent  of  i      and  go  the  right  side  of  (lo5)  is  less 
than  (lo2);  this  contradicts  our  assertion  that  (1.32)  implies 
(1.35).   The  resolution  of  this  contradiction  is  that  the 
solution  does  not  exist  for  values  of  T  violating  ( 1-35 ) • 
This  completes  the  proof  of  Theorem  1.1. 
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2 .   Approximate  conservation  laws  and  approximate  characteristics^ 
We  start  by  outlining  briefly  the  theory  of  shock  waves,  rare- 
faction waves  and  their  interactions.   For  details  see   [  ^]   and 
[2]. 

A  discontinuous  function  which  is  a  solution  in  the  sense  of 
distributions  of  the  conservation  laws  (l.l)  satisfies  along  the 
curve  of  discontinuity  the  jump  relations 

(2.1)  s[u]  =  [f],  s[v]  =  [g] 

where   [u]   denotes  the  difference  of  the  value  of  u   on  the  right 
and  on  the  left  side  of  the  discontinuity,  and   s   denotes  the  speed 
with  which  the  discontinuity  propagates.   It  is  well  known  that 
physically  meaningful  discontinuous  solutions  satisfy  in  addition 
the  entropy  condition;  under  the  normalization  (1.4)  this  takes  the 
form 

(2.2)^  [w]  <  0 

for  discontinuities  propagating  to  the  right   (s   positive), 

(2.2)_g  [z]  <  0 

for  discontinuities  propagating  to  the  left.   Discontinuous 
solutions  satisfying  (2.1)  and  (2.2)  are  called  shocks. 

It  has  been  shown  in  [ 4  ]   that  for  weak  shocks  propagating 
to  the  right   (0  <  s) 
(2o3)  [z]  =  0[w]^ 

while  for  weak  shocks  propagating  to  the  left   (s  <  O) 

[w]  =  0[z]^. 
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For  the  equations  of  gas  dynamics  both   z   and  w  decrease 
across  shocks  propagating  in  either  direction: 

(2.4)         [z]  <  0,   [w]  <  0  across  shocks. 

We  stipulate  that  this  condition  be  fulfilled  for  the  systems 
studied  in  this  paper. 

It  follows  from  (2.1)  and  (2.2)  that  for  weak  shocks 

([z],  [w]  small)   the  characteristic  speed  in  the  direction  in 

which  the  shock  is  moving  is  greater  than  the  shock  speed  l3ehlnd 
the  shock  and  less  in  front  of  it.   I.e., 

p   <  s  <  p„    when    0  <  s 

(2.5) 

Aj,  <  s  <  A.    when   s  <  0; 

here  the  subscript  I      and  r   refer  to  the  state  to  the  left 
and  to  the  right  of  the  shock. 

The  simplest  and  most  basic  initial  value  problem  is 
the  Riemann  problem  where  the  initial  value  has  one  constant 
value   a   for   x  negative  and  another  constant  value  b   for 
X  positive.   It  is  not  hard  to  show  that,  at  least  when  a 
and  b   are  close  enough,  this  initial  value  problem  has  a 
solution,  depending  only  on  x/t;    this  solution  may 
contain  a  shock  propagating  in  either  or  both  directions. 
The  solution  is  constant  in  certain  sectors  and  in  those 
sectors  where  it  isn't, either   z   or  w  is  constant.   These 
regions  are  called  centered  rarefaction  waves;  they  propagate 
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either  to  the  right  or  to  the  left,  and  their  boundaries  are 
characteristic  lines.   Figure  2.1  presents  one  of  the  four 
possibilities: 


left 
rarefaction  wave 


right  shock 


Fig.  2.1 


Since  the  right  edge  of  a  centered  rarefaction  wave 
moves  faster  than  the  left  edge,   [p]   and  [a]   are  positive 
for  right,  respectively  left,  rarefaction  waves, where   [   ] 
denotes  the  difference  of  the  value  of  the  function  in 
brackets  on  the  right  edge  of  the  rarefaction  wave  and  its 
value  on  the  left  edge.   In  view  of  the  normalization  given 
in  (1.4)  these  relations  are  equivalent  with 


(2.6) 


0  <  [w],   0  <  [  z] 


respectively. 

We  indicate  on  figure  2.2  the  relative  positions  in 
the  w,   z  plane  of  the  constant  states   a,   m  and  b 
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occurring  in  the  interaction  pictured  in  figure  2.1; 


/N 


^ 


w 

Fig.  2.2 

Relations  (2.2)  ,  (2-3).  and  (2.4)    for  the  left  shock 

connecting   a  and  m  express  the  fact  that   m  lies  below 
a  and  slightly  to  the  left  of  a,   while  the  constancy  of   z 
in  a  rarefaction  wave  moving  to  the  right  and  relation  (2.6) 
express  the  fact  that   b   and  m  lie  on  the  same  horizonal 
line,   b   to  the  right  of  m. 

Figure  2.3  shows  two  shock  waves  moving  to  the  right 
which  coalesce  to  form  another  shock  wave  moving  to  the  right 
and  a  rarefaction  wave  moving  to  the  left: 


A 


t 


-> 


X 


Fig.  2.3 


The  explanation  of  this  phenomenon  lies  in  relations  (2.2), 
(2.3)  and  (2.4),  as  illustrated  by  figure  2.4: 
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A 


■>■ 


w 


Fig.  2.4 


We  shall  denote  by   [z]  and  [w]   the  magnitude  of  a 
wave  propagating  to  the  right  and  left  respectively;  according 
to  (2.2)  and  (2.6)  the  magnitude  of  a  shock  wave  Is  negative, 
of  a  rarefaction  wave  positive. 

We  shall  denote  by   (a,b)    and   (a,b)    the  magnitude 
of  the  left  and  right  waves  which  arise  In  the  resolution  of 
the  Rlemann  problem  with  state   a  on  the  negative  axis,   state   b 
on  the  positive  axis.   We  shall  say  that   a   can  be  joined  to  b 
through  the  waves   (a,b),   and   (a,b)  . 

Suppose  that  a  can  be  joined  to  b  through  a  single 
rarefaction  wave  moving  to  the  left;   let   c   be  another  state 
which  can  be  joined  to  b   through  a  rarefaction  wave  moving 
to  the  left  and  some  other  wave  moving  to  the  right.   Then 
clearly  a  can  be  joined  to   c   through  a  left  rarefaction 
wave  which  is  the  sum  of  those  joining   a  to  b   and  b   to 
c,   and  the  right  wave   (b,c)  . 

An  Important  generalization  of  this  simple  observation  la 
given  In  a  the  following  theorem  of   [ 2  ] : 
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Theorem:   Let   a,  b,  c   be  any  three  nearby  states; 


then 


(2.7)^    (a,c)^  =    (a,b)^  +  (b,c)^  +  error 

and 

(2.7)^    (a,c)^  =  (a,b)^  +  (b,c)^  +  error. 


The  error  is  of  the  form 

(2.8)  error  =  0(Ia-bJ  +  |b-c | )Q(a,b, c) 

where   Q  is  a  measure  of  the  interaction  of  waves  connecting 
a  t_o  b  with  those  connecting  b   to_  c .   The  precise 
definition  of  this  notion  is  as  follows: 

A  wave  connecting   a  to  b   interacts  with  those 
waves  connecting  b   to   c   which  it  would  intersect.   Thus 
the  right  wave  connecting   a  to  b   interacts  with  the  left 
wave  connecting  b   to   c   while  the  left  wave  connecting 
a  to  b   does  not  interact  with  the  right  wave  connecting 
b   to   c.   On  the  other  hand  the  two  left  waves  interact  with 
each  other  if  and  only  if  one  of  them  is  a  shock  wave,  and 
similarly  for  the  right  waves. 

The  numerical  measure   Q  is  defined  as  follows: 

(2.9)  Q(a,b,c)  =  ^  |(a,b)^(b,c)jl, 

the  summation  being  over  all  pairs  of  interacting  waves. 

We  describe  now  briefly  the  difference  scheme  used 
in   [2  ]   to  construct  approximate  solutions.   An  approximate 
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solution  Is  exact  in  each  strip  ns  j<  t  <  (n+l)s;   at 
t  =  ns   the  solution  Is  constant  in  each  interval 
(in-2)r  <  X  <  mr.   The  ratio  r/s   is  chosen  to  be  greater  than 
Max(pjA|);  this  guarantees  that  the  waves  generated  at  the 
points  of  discontinuity  x  =  (m-2)r   and  m  r   don't  interact 
with  each  other  in  the  above  strip,  so  that  the  exact  solution 
there  consists  of  constant  states  separated  by  centered 
rarefaction  waves  or  shocks,  see  figure  2.5".: 
n+1^ 


n 


m-2 


m 


m+2 


m+4 


Fig.  2.5 


We  shall  call  these  waves  elementary  waves  of  the 
approximate  solution  under  discussion. 

At   t  =  (n+l)s   the  exact  solution  is  not  piecewise 
constant;  an  approximate  piecewise  constant  solution  is 
obtained  according  to  the  following  recipe: 

Pick  a  random  number  '^]^^]_>   equidistribution 
in   (-1,1),   and  set  the  value  of  u  at  time   (n+l)s   in  the 
interval   (m-l)r  <  x  <  (m+l)r   to  be  equal  to  the  value  of 
the  exact  solutions  constructed  in  the  strip  at   t  =  (n+l)s 
and  X  =  (^+°'n+l^''"* 

Equidistribution  means  that  for  every  subinterval   I   of   (-1,1) 
lim  ^^[|^^^  =  -I  length  of   I, 
where   S(N,l)   is  the  number  of  a  ,   n  <  N  which  lie  in  I. 
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m-2-fa  m+a  ^  m+2+a  , 

n+1  n+1  n+1 


n+1  x=^ 5t-7 r e:^&=7 >^ 5^ 


n 


m-2  m  m+2 

Fig.  2.6 

Iri   [  2  ]   the  following  theorem  was  proved: 

Suppose  the  initial  values  satisfy  the  following 
conditions : 

i)   The  oscillation  e   of  the  initial  data  is  small  • 
ii)   T)  =  e   (total  variation  of  data)    is  small. 

Then  there  exists  a  sequence   r   and   s   tending  to 
zero,   r/s   constant  subject  to  the  above  inequality, 

'  -■■■  _  ■  ■      ■      ■      ■  ■  ■    ■  ^y  ■  ■  ■  ■  - 

such  that  the  approximate  solutions  tend   to  a  weak  solution 
of  the  conF^ervation  equations  (l.l)  for  almost  all  random 
sequences   ^  i   .   The  oscillation  of  this  solution  is   0(e) . 

A  basic  tool  in  the  proof  of  this  theorem  is  a  careful 
analysis  of  the  interaction  of  elementary  waves.   This  interaction 
takes  place  in  the  diamond -shaped  region  A     whose  vertices 
S,  W,  E,  N  are 

S  =  {"i4^n-l  '^"^^  '      ^  =  (m-HX^^^,n+l) 
W  =  (m-l+a^,n)   ,     E  =  (m+l-KX^^^n)  . 

In   [ 2  ]   the  random  sequence  depended  on  m  as  well  as   n; 
it  is  easy  to  show  that  one  can  dispense  with  the  randomization 
in  m. 
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The  point   (m,n)   is  called  th-e  center  of  A 

ill  }  l  i- 

The  collection  of  all  diamonds  covers  the   x,t   plane,   t  >  0. 


N 


n+1 


n  W   x<    o       A  _   ^x  E 


n-1 


m-1 


Fig.  2.7 


The  elementary  waves  originating  at  time  n-1  at  the 

points   m-1   and  m+1   enter  A    ,   interact  there,  and  give 
^  m,n^  '  ■= 

rise  to  two  waves  emerging  from  the  center   (m,n)   and 

propagating  to  the  right  and  left,  respectively. 

According  to  (2-7)  the  magnitude  of  each  of  the  waves 

emerging  from  A   equals  the  sum  of  the  magnitudes  of  the  waves 

belonging  to  the  same  family  which  enter  A,   plus  an  error 

of  size  (2.8).   However  when  two  waves  of  the  same  family  enter 

A  and  one  of  them  is  a  rarefaction  wave,  the  other  a  shock, 

as  illustrated  in  figure  2.7,  cancellation  takes  place.   The 

amount   C  .  ,  ,   of  right  waves  cancelled  in  A   is  defined  to 
right       ^ 

be 

(2.10)  ^right^^^  =  I  (|a|+|p|-|a+p|) 
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where  a   and   P   are  the  magnitudes  of  the  right  waves  which 
enter   A. 

We  can  summarize  these  facts  as  a 

Conservation  law  for  diamonds; 

Denote  by  E— .  u+(A)   the  sum  of  the  magnitudes  of  all 
"^        right  ^  ° 

rarefaction,  respectively  shock  waves  moving  to  the  right 
which  enter  A,   and  denote  by  L— .  j„^,  (A)   the  magnitude   of 
the  right  wave  which  leaves   A.   Then 

(2.11)     Li.  K+(A)  =  Ei.  ,  +(A)  T  C  .  v.+  (A)  +  0(£)Q(A) 
\     'r    right  ^      right  ^  '  +  right  ^     y^\    > -^x    i 

where   C  .  ,  , (A)   denotes  the  amount  of  right  wave   cancelled 
right  ^ 

in  A,   e  is  the  oscillation  of  the  solution  in  A   and  Q(A), 

defined  by  (2.9),  measures  the  interaction  of  all  waves  entering 

A.  A  similar  law  holds  for  left  waves. 

Let  J\,   be  the  union  of  a  collection  of  diamonds, 

E— .  u,j.(A  )   and   L— .  ,  ,  (^)   the  sum  of  the  magnitude  of  the 
right^-*'*-'         right^^"*"  ^ 

right  waves  entering  and  leaving  y\.  ,   respectively.   By 
adding  up  the  conservation  laws  (2.11)   for  the  diamonds 
constituting  J^     we  get  the  following 
Conservation  Law: 

Where  C  .  ,  ,  (A)   is  the  sum  over  all  diamonds  of  ^  of 
right  ^"^  *• 

the  amounts  of  right  waves  cancelled  in  each  diamond  A;   e   is 
an  upper  bound  for  the  oscillation  in  each  diamond,  and 
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(2.13) 


Q(A)  =  XL       Q(^) 
A  in  A 


A  similar  conservation  law  holds  for  left  waves. 

Let   I   be  any  interval  on  the  initial  line, 
the  domain  of  determinacy  of   I,   see  figure  2.8: 


A(I) 


V. 


Fig.  2.8 


Denote  by  V(l)   the  total  variation  of  w  and   z  along  I> 
it  has  been  shown  in  [2  ] ,      inequality  on  top  of  page  709, 

that  if 

eV(  I )  <_  const, 

then 


(2.14) 


Q(/\(I))  <  2V(I)^. 


This  estimate  is  basic  for  this  paper. 

Since  rarefaction  waves  have  positive  magnitude. 


+ 


+ 


L   are  nonnegative.   Thus  we  deduce  from  (2.12)    that 


(2.15) 


C^(/\(I))  <  E^{A{1))    +   0(e)Q(A(l)) 


1 


Note  that  the  Q   in   [  2  ]   is  different  from 
paper,  which  is  denoted  there  as   D. 
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in  this 


The  waves  entering   y\.(l)   are  just  the  elementary  waves 
emanating  from  the  initial  interval   I,   and  therefore  we  can 
estimate  the  sum  of  their  magnitudes  in  terms  of  the  total 
variation  of  w  and  z  along  I: 

(2.16)  E+.g^^(A(l))  +  E+^^^(/\(l))  <  0(V(I)) 

Using  (2.16)  and  (2.l4)  to  estimate  the  right  side  of  (2.I5) 
we  get 

(2.17)  Cl^ight^'^^^^^  +  ^left(^(^^)  -  0(V(I))  +  0(ev2(i}). 

In  addition  to  the  conservation  law  for  domains 
composed  of  diamonds  we  need  one  for  domains  boiinded  partly  by 
characteristics  and  shocks.   For  this  purpose  we  define  for 
each  approximate  solution  two  families  of  approximate  characteristic 
curves,  propagating  to  the  left  and  right  respectively,  as  follows: 

In  each  strip   (n-l)s  <  t  <  ns,   where  the  approximate 
solution  is  exact,  an  approximate  characteristic  is  a  straight 
line  segment  "^ ^  ~    '      which  is  an  exact  characteristic  or 
shock  and  which  originates  at  time   t  =  (n-l)s   from  the 
center  of  some  diamond.   The  continuation  of  this  characteristic 
is  a  similar  line  segment  "X.  located  in  the  next  strip  and 

issuing  from  the  center  of  that  diamond  A   in  which  the  upper 
endpoint  of  y^^^'^^      lies: 
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n+1 


n 


n-1 

Fig.  2.9 

The  slope  of  %  has  to  be  determined  so  that,  except 

for  cancellations,  there  is  a  conservation  of  waves  on  both  sides 
of  the  characteristic.   We  will  make  this  more  precise  after  giving 
the  rule  in  each  of  the  six  possible  situations  depending  on  the 
nature  and  relative  strength  of  the  two  waves  of  the  same  family 
entering  A;   the  trivial  case  when  only  one  wave  of  the  family 
under  discussion  enters  A  is  subsumed  under  our  classification 
if  we  regard  the  second  wave  as  having  strength  zero. 

In  figure  2.10  only  right  waves  are  shown,  two  entering, 
one  leaving. 

Since  each  wave  could  be  a  shock  or  a  rarefaction,  there 
are  eight  possible  combinations.   According  to  (2.7)  the  magnitude 
of  the  wave  leaving  is  the  sum  of  the  magnitudes  of  the  waves 
entering,  plus  an  error  term;  this  shows  that  if  we  disregard  the 
error  term,  then  if  the  two  entering  waves  are  both  shock  or  both 
rarefaction,  the  wave  leaving  is  also  a  shock  or  rarefaction. 
Thus  out  of  the  eight  possibilities  only  six  are  dynamically 
permissible. 

'yC      ~'^        and   7f     are  indicated  by  dotted  lines;  for  the 
sake  of  definiteness  we  assume  that  y^^  ~        enters  from  the  left. 
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i) 


lii) 


vi) 


In  cases   i)  -  iii)   X 


(n) 


is  the  shock  wave  leaving 
v(n) 


A.   In  cases   iv)  -  vi)   the  choice  of  "pC"^'      ^^   dictated  by 
the  following  principle: 

Rarefaction  waves  cannot  cross  approximate  characteristics. 

In  case  iv)  this  means  that  the  magnitude  of  the  rarefaction 
wave  leaving  A   to  the  right  (left)  of  X  equals  the  sum  of 

the  magnitudes  of  the  rarefaction  waves  entering  A  to  the  right 
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(left)  of  7l^^~^\ 

( n) 
In  case  v)  -/^j  is  the  left  edge  of  the  rarefaction 

wave  leaving  A. 

In  case  vi )  there  are  two  subcases:  if  the  magnitude  of 

the  rarefaction  wave  leaving  A   is  greater  than  the  magnitude 

of  the  wave  which  enters   A   to  the  left  of   7i     >  we  choose 

A.     so  that  the  magnitude  of  the  wave  leaving  A   to  the 

(n) 
left  of   7v     equals  the  magnitude  of  the  wave  entering  A 

to  the  left  of  7^^  ~   .   Otherwise  we  choose   Tt,     to  be  the 

right  edge  of  the  rarefaction  wave  leaving  A. 

According  to  (2.7)  the  magnitude  of  the  rarefaction  wave 

leaving  A   differs  by   0(e)Q(A)   from  the  sum  of  the  magnitudes 

of  the  waves  entering  A;  it  can  happen  therefore  in  case  iv) 

and  in  the  first  subcase  of  vi )  that  X     prescribed  above 


lies  outside  the  rarefaction  fan  leaving  A;   in  such  a  case  we 

.(n-1) 


redefine  /V     to  be  the  right  edge  of  the  wave 


Suppose  two  different  approximate  characteristics  X 
and   ■^^^~^'      enter  the  same  A,  Tp^^''^'      to  the  right  of  -^^ 
It  is  not  hard  to  see  that  our  prescription  puts  7(^     to  the 
right  of  TV   ;   thus  approximate  characteristics  of  the  same 
family  cannot  cross.   Of  course  it  can  happen  -  in  fact  one  of 
the  main  conclusions  of  this  paper  is  that  it  must  happen  - 
that  ^^     coincides  with  7t   .   Thus  initially  distinct 
approximate  characteristics  may  coalesce . 

Connect  the  endpoint  of  X       to  the  beginning  of 
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/{        ;      the  resulting  plecewise  linear  curve  divides   A   intc 
two  halves: 


left  a/2 


right  a/2 


Fig.  2.11 

In  each  half,  denoted  as   A/2,  the  following  conservation  law 
for  right  waves  holds: 

(2.18)^    L"^(A/2)  =  E^(A/2)  -  C^(A/2)  +  0(e)Q(A/2) 

and 

(2.18)"     S{A/2)  +  L"(A/2)  =  E"(A/2)  +  C"(A/2)  +  0(e)Q(A/2) 

where : 

Li(A/2)   and  &i(A/2)   denote  the  amount  of   + 
right  wave  entering  and  leaving  A/2   through  that  part  of  its 

boundary  which  is  part  of  the  boundary  of  A.   If  tC 

(n) 
or  X     is  a  shock  wave,  we  count  it  as  belonging  to  neither 

half  of  A;   for  this  reason  we  introduce  the  nonpositive  quantity 

S{A/2)   to  denote  the  amount  of  shock  wave  which  leaves  A/2   and 

becomes  part  of  the  interior  boundary. 
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C— (a/2)   are  nonnegatlve  quantities  and  denote  the  absolute 
amount  of  shock  and  rarefaction  waves  respectively  cancelled  in 
a/2.   These  are  no  longer  equal;  nevertheless  the  following  is 
true: 
(2.19)      Q( right  A/2)  +  Q( lef t  A/2)  -  Q(A) 

(2.19)^     C^( right  A/2)  +  C^( lef t  A/2)  =  C(A) 
and 

(2.19)"      C"(right  A/2)  +  C~( left  A/2)  ^  C(A). 

where   right  A/2   and   left  A/2   denote  the  two  halves  of  A. 
The  reason  why  the  sign  of  inequality  may  hold  in  (2.19)~   is 
that  some  of  the  cancellation  may  involve  a  shock  on  the 
characteristic  and  that  is  not  counted  as  part  of  either  right 
A/2   or   left  A/2 . 

The  conservation  law  may  he  verified  in  each  half  of  A 
in  all  six  cases. 

If  there  are   2,  3,  ...,  m  approximate  characteristics 

crossing  A,   they  divide   A   into  m  +  1   parts;  in  each  part 

a  conservation  law  like  ( 2 . l8 ) ,  (2.19)  holds.   Let  Q.   be  the 

J 
th 
amount  of  interaction  present  in  the  j    fraction,   j  =  1,  ..., 

m  +  1.   It  can  be  verified  that 

XI  Qj  =  0(Q(A)). 

Let  /\    be  the  union  of  a  collection  of  diamonds  and 
diamonds  cut  by,  say,  right  approximate  characteristics. 
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Adding  up  the  conservation  laws  (2.11)  for  the  diamonds  and 
(2.18)—  for  fractional  diamonds  we  get  the  following 
conservation  law: 

(2.20)+     L'^(A)  =  E+(A)  -  C  +  (A)  +  0(e)Q(A) 
and 

(2.20)"     L"(A)  +  S(^A)  =  E"(A)  +  C"(A)  +  0(e)Q(A) 

where   S(Syl)   denotes  the  amount  of  shock  wave  which  leaves 
J\.    and  "becomes  part  of  its  boundary. 

We  shall  use  these  conservation  laws  to  derive  an  estimate 
for  the  total  variation  of  w  and   z   along  an  approximate 
characteristic   7(  .   We  estimate  first  the  total  amount  of  shock 
entering   ^  and  then  the  total  amount  of  left  wave  crossing   ^. 


t 

o 


Fis.  2.12 
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Denote  by  p   the  point  where  the  characteristic  segment  at 
the  upper  end  of  'X   intersects  the  boundary  of  the  diamond 
from  which  it  is  emerging.   Construct  the  curve   C  which 
connects  the  point   p   to  a  point   q   on  the  line   t  =  t 
on  which   7l  starts  so  that   <r  consists  of  NE   sides  of 
diamonds  along  which   t   decreases  from  p   to  q  ,   see 
figure  2.12.    Denote  by  V\  the  union  of  diamonds  and  half 
diamonds  between  X  ^•i^d   C  and  lying  above   t  =  t  .   The 
only  waves  which  enter  _/V  enter  through  the  part  of  the 
boundary  labeled  <y         which  lies  on   t  =  t  ;   therefore  we 
conclude  from  the  conservation  law  (2.20)  that 

(2.21)  |S(X)1  1  |e;(o'^)|  +  0(e)Q(A). 

Denote  by  J\.        the  union  of  whole  diamonds  which  are 
totally  contained  in  /\.   and  apply  the  conservation  laws  (2..  12) 
for  left  waves  to  ^V  •   The  left  waves  which  leave  W   on 
the  boundary  adjacent  to  X  are  the  ones  which  enter  the 
diamonds  in  which   %  is  situated.   We  call  their  s\am  the 
total  amount  of  left  wave  crossing  X^      ^^^   denote  it  by 
X— (X  )  •    Since  all  the  left  waves  entering   /\   do  so  through 
cr   we  conclude  from  the  conservation  law  (2.12)   that 

(2.22)  X"'(X)  +  |X-(X)|  lE+gft^'^o^  +  l^left^'^o^l  +  0(e)Q(A^) 

Next  we  make  use  of  the  last  two  estimates  to  estimate 
the  total  variation  of  w  and   z   along  either  side  of   X- 
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It  is  easily  verified  by  inspection  of  all  cases  described 
in  figure  2.10  that  the  change  in  w  across  right  waves 
between  the  points  where  tC   enters  and  leaves   A   does  not 
exceed   S(A/2)  +  C^(A/2)  +  0(e)Q(A/2).   A  further  change  in 


w 


is  due  to  the  left  shock  waves  crossing  X  ?  but  the 


2 
magnitude  of  this  is   0(e)   times  the  strength  of  that  shock 

wave.   Since  we  have  already  obtained  bounds  for  the  total 

amount  of  shock  entering  ~/i,      the  total  amount  of  left  wave 

crossing   7C^   ^-^^  '^^^   total  amount  of  cancellation  in  the 

diamonds  containing   TC  j   'we  obtain  an  estimate  for  the  total 

variation  of  w   on  either  side  of   7^..   A  similar  estimate 

holds  for  the  total  variation  of   z   along  either  side  of   7i • 

As  an  easy  corollary  we  estimate  the  total  variation 

of  the  speed  of  the  segments   7t     which  make  up   tC*   If 

X^^   lies  in  a  rarefaction  wave,  its  speed  '^^  equals 

p(w, z),   w  and   z   being  the  value  of  the  variables  at  the 

( n )  (  n) 

point  where   tC     leaves   A;   if   7^^     is  a  shock,  its 

speed,  according  to   [  '^  ]  ^   is  approximately  the  average  value 

of   p   on  the  two  sides;  more  precisely 

.  i^\         p(w  ,z  )+p(w  ,z  )  p  p 

(2.25)   7C^    =  = = ^-^  +  0((w  -w  r    +  (z  -z  )'^). 

2  +  + 

Of  course  (2.25)  is  valid  in  both  cases.   It  follows  easily 

.  ( n) 
from  (2.25)  that  the  total  variation  of   7^     along  -v     is_ 

0( total  variation  of  w  and   z   along  both  sides  of  X^* 

We  need  one  more  conservation  law,  and  for  this  purpose 

we  introduce  the  concept  of  an  approximate  time-like  curve. 
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We  recall  from  the  theory  of  hyperbolic  equations  that  a  curve 
X  =  x(t)   is  called  time-like  if  its  speed  of  propagation  x 
lies  between  the  two  characteristic  speeds: 

(2.24)  A  <  X  <  p. 

We  assume  inequality  (2.24)  to  hold  for  all  the  left  and  right 

speeds   A   and   p   at  any  point  in  the  flow  field. 

Given  an  approximate  solution  we  construct  an  approximation 

to  a  given  time-like  curve  as  follows:   t      is  a  straight  line 
^  appr  ^ 

segment  in  each  strip   (n-l)s  <  t  <  ns,   whose  speed  equals   x(ns), 

starting  at  the  center  of  that  diamond  which  contains  the  endpoint 

of  the  previous  segment. 

Let   T   be  an  approximate  time-like  curve;  denote  by  a(x) 

the  union  of  all  diamonds  which  contain  points  of  t.   It  is  an 

Immediate  consequence  of  (2.24)  that  left  waves  enter  a(t)  only 

from  the  right  and  leave  it  only  on  the  left;  the  converse  is 

true  for  right  waves.   The  amo-unt  of  right  and  left  waves  entering 

A(t)   is  defined  to  be  the  amount  of  right  and  left  wave  crossing 

t;   these  quantities  are  denoted  by  X—.  ,  ,  (x)   and  X^  ^^.(x) 
'        ^  "^    right  left 

respectively,  where  the  superscript  +  indicates  rarefaction 

and  shock  waves  . 

The  conservation  law  in  question  concerns  domains  bounded 

by  an  approximate,  say,  right  characteristic   X,   an  approximate 

time-like  curve   x  and  an  interval  c       on  a  line  t  =  const, 
o ' 

It  follows  from  the  conservation  law  (2.12)  that  the  amount  of 
wave  which  actually  crosses   x,   i.e.  emerge  on  the  other  side 
of  A(x)   is  X(x)   minus  the  amount  of  waves  cancelled  in  A(x), 
modulo  0(e)Q(A(x)). 
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as  indicated  schematically  in  figure  2.13: 


o 


Fig.  2.13 

Consider  the  region  /\.   consisting  of  all  diamonds  which  are 
entirely  contained  between  yL,      i     and  c  .   According  to  the 
foregoing,  left  waves  enter  y\.    only  through  the  bottom  and 
right  boundaries,  and  leave  it  only  through  the  left  bo\indary. 
Applying  the  conservation  law  (2.12)  we  conclude: 


(2.25) 


l^left^^o^ 


-  '^left^^^l  ^  ^^^^    ^   0(e)Q(A), 


where 


X 


is  an  abbreviation  for  X-,  ^.  +  |X-,  ^, 


left 
¥e  could  define  approximate  space-like  curves   o- 


analogously  but  that  is  not  necessary;  we  define  instead  an 

approximate  space-like  curve  to  be  a  continuous  curve  which 

is  a  union  of  WN   or  WS   sides  of  diamonds.    The  amount  of 

right  and  left  wave  crossing   o-  is  defined  in  the  obvious 

fashion  and  is  denoted  as  X-i  ^,  (cr)   and  X-i.  ,  .(C).   We  abbreviate 

left  right 


38 


by  the  symbol   |x(o')|   the  total  amount  of  wave  crossing   C: 

\x{<y)\    -X+^^^(a)  +  |x-^^^(o')|  +X+    (a)  +  |x;,g^^(cr) 


The  total  amount  of  wave  crossing  a  space-like  curve 
o'  is  almost,  but  not  quite,  equal  to  the  total  increasing  and 
decreasing  variations  of  w  and   z   along   o";   the  discrepancy 
is  due  to  the  fact  that   z   changes  across  a  right  shock  wave, 
although  only  by  an  amount  proportional  to  the  cube  of  the  shock 
strength.   From  this  estimate  for  the  discrepancy  we  deduce  the 
following  relations: 

Increasing  variation  of  w  along  cr  = 

(2.26)         =  X  +  .g^^(or)  +  O(e^)|x(0')|. 

Decreasing  variation  of  w  along  a   = 

=  ^right^^^  +  O(e^)|x(or)|. 

Similar  relations  hold  between  the  variation  of  z  and  the  total 
amount  of  left  wave  crossing  cr. 

We  summarize  for  easy  reference  the  various  estimates, 
see  (2.17),  (2.21),  (2.22),  derived  with  the  aid  of  our  conser- 
vation laws: 

Theorem  2.1:   The  quantities  listed  below  are  bounded 
by   0{Y{<y))    +   0(e)V(o')  ,   where  V(a)   is  the  total  variation 
of  w  and   z   along  a  space-like  curve  cr     whose  domain  of 
determinacy  includes  the  pointset  in  question; 
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(a)   C  .  ,  j_iA)    +   C-,  ^.iA),    the  amount  of  wave  cancelled 
^      right       left     

in  a  union  A   of  diamonds. 

( B )  S ( X) ,   the  amount  of  shock  wave  which  enters  a 
characteristic . 

(C)  The  total  variation  of  w  and   z   along  either 
side  of  a  characteristic   X- 

(D)  The  total  variation  of  the  speed  ^  of  a 
characteristic . 

(E)  Q(A),   the  total  amount  of  interaction  in  A. 
¥e  derive  now  some  further  estimates  on  the  strength 

of  shocks: 

The  analysis  which  led  to  the  conservation  law  (2.l8) 
shows  that,  mod  0(e)Q(A),   the  strength  of  a  shock  crossing 
A  can  Increase  at  most  by  the  amount  of  shock  which  enters 
A,   and  can  decrease  at  most  hy  the  amount  of  rarefaction  waves 
cancelled  in  A. 

Let   X  be  an  approximate  characteristic  and  denote  by 
y\.      the  union  of  all  diamonds  through  which  X   cuts  during  the 
time  t.  :!  t  ^  t„.   Denote  by  Str'?C(t)   the  strength,  i.e.  the 
absolute  value  of  the  magnitude,  of  the  shock  on  %    at  time   t. 
The  estimate  given  above  on  the  increase  or  decrease  of  the 
strength  of  X   'u-pon  crossing  each  diamond  yields  this  inequality; 

(2.27)       StrX(t^)  -  C(A)  +  0(e)Q(A)  1  StrX(t^)  < 

<  StrX(t^)  +  |e~CA)|  +  0(£)Q(A), 
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E  CA)   being  the  amount  of  shock  entering 

Using  the  definition  (2.9)  of  interaction  Q  we  can 
establish  the  following  lower  bound  for  Q(/\)  - 

(2.28)  Str^.^|E-(A)|  1  Q(/A) 

where   Str  .    is  the  minimum  strength  of  X  in   (t.,t„). 

We  give  now  an  estimate  of  the  amount  of  shock  E~(A) 
which  enters  /( .      Let   I   be  an  interval  at  time   t  =  t. 
containing    ;^(t.);   suppose  that  there  are  approximate  right 
characteristics  Xp      ^^^   X    issuing  from  the  endpoints  of 
I   which  do  not  run  into  X   fo^  ^  1  '^f      Denote  by  y\        the 
set  of  points  between  Xn      s-'^d  X    ,      t.  <  t  <  t„.   According 
to  the  conservation  law  (2.20)~ 

(2.29)  lL"(Ax)|  1  |E"(/lj)|  +  0(e)Q(Aj). 

Since  j\-^      is  bounded  by  approximate  right  characteristics, 
shock  waves  enter  it  only  through   I,   i.e. 

E'(/^j)  =  X^(I). 

On  the  other  hand  since    Xp      and  X        do  not  run  into  X, 
X(t^)   is  one  of  the  waves  leaving  yV;   so  we  have  from  (2.29) 
that 

(2.30)  StrX(t^)  1  |X^(I)|  +  0(e)Q(Aj). 
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3 .    The  construction  of  exact  solutions  as  a  limit 
of  approximate  solutions. 

The  basic  estimate  of   [  2  ]  is  that  the  approximate 
solutions   u,   described  in  Section  2   satisfy  the  inequality 

(3.1)        Var  u,  (x,t)  <_  const  Var  u,  (x,0). 

Furthermore,  as  was  shown  in  [  2  ] ,  this  inequality  implies 
that  the  approximate  solutions   u,  (-jt)   are  uniformly 
continuous  in  the  L-,°^   topology.   It  follows  from  this, 
with  the  aid  of  Kelly's  theorem,  that  a  suitably  selected 
subsequence  of   1  ^v  }  converges  in  the  L-,     topology 
with  respect  to   x,   to  some  limit   u,   uniformly  in   t. 
It  was  shown  in   [  2  ]   that  for  almost  all  choice   of 
the  sequence  of  random  points  |^ct„  3  '^^^   limit   u   is  a 
weak  solution  of  the  originally  given  conservation  laws. 

The  limit   u   is  of  bounded  variation  in  x   for 
each  t;  it  follows  that  for  each  t,   u(x,t)   is  a 
continuous  function  of  x   except  for  a  denumerable  set 
of   X.   It  is  not  hard  to  show,  using  Fubini's  theorem, 
that  the  set  of  discontinuities  has  measure  zero  in  x,t. 

The  purpose  of  this  section  is 

a)  To  construct  characteristics  and  shocks  for 
exact  solutions. 

b)  To  show  that  the  slope  of  such  a  characteristic 
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equals  the  characteristic  speed  at  that  point,  and  that  the 
speed  of  a  shock  satisfies  the  Rankine-Hugoniot  equations. 

c)  To  prove  approximate  conservation  of  variation  for 
regions  bounded  by  space-like,  time-like  and  characteristic 
arcs . 

We  shall  construct  characteristic  curves  yC    for  an 
exact  solution  u  as  limits  of  approximate  characteristics 
X,       constructed  in  Section  2  for  the  approximate  solutions 
u,   of  which  u   is  the  limit: 

Lemma  3»1:   Let   x  =  /^,  (t)   be  an  approximate  character- 
istic, whose  initial  point  {%..   (t,  ),t,  )   tends  to   (x  ,t  ). 

K   K     K  O   O 

Then  we  can  select  a  subsequence  of  k   such  that  ^|^(t) — ^;^(t) 
uniformly  on  every  bounded  time  interval.   The  limit  %[\^)      is 
Lipschitz  continuous. 

Proof;   By  construction,  %,       satisfies 


(3.2)     \X^{t^\-X^{t^)\    <    (V^k^/lH-^2l  +  2^ 


k 


Since  in  the  difference  scheme   ^v/^i,   is  independent  of  k, 
(3.2)  gives  a  uniform  bound  on  the  oscillation  of  /^(t); 
the  existence  of  a  uniformly  convergent  subsequence  follows 
then  as  in  the  argument  for  the  Arzela-Ascoli  theorem. 

(3.2)  implies  that   ^v/^i^^   is  a  Lipschitz  constant  for 
the  limit  function   X. 

We  recall  from  Section  2  that  an  approximate  character- 
istic consists  of  straight  line  segments,  and  that  the  point 

th 
of  origin  of  the  n   segment  depends  on  the  choice  of  a 

random  parameter  a  .   The  following  lemma  is  crucial  for 
^  n  ^ 
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showing  that  for  a. a.  choices  of  the  set  of  random  parameters, 
the  limit  characteristic  curve   has  the  correct  speed  of 
propagation. 

Lemma  3-2:   a)   Let  the  exact  characteristic   X  be 
the  limit  of  the  approximate  ones  '^y'      ^^^    suppose  that 
for  all   k   large  enough  and  all   t   in   [t^,t*]   the 
speeds  X.y      satisfy 

(3.3)^  a  <  i^  I  b- 

Suppose  that  the  sequence  \^-\.\     of  random  numbers  is 
equidistributed;  then  the  speed  of  the  limit  characteristic 
also  lies  between  a  and   b: 

(3.3)  a  <  i  <  b 

b)   Suppose  that  %.    and   ^   are  two  characteristics, 
limits  of  the  sequences  %  ;,   3-nd  ^,   respectively. 
Suppose  that  for  all   k   large  enough  and  all   t   in 
[t^,t*]   the  speeds  of  the  approximate  characteristics 
satisfy 

(3.4)  "^k  -  ^k  +  "^  • 

Suppose  that  the  sequence   ^a  f   is  equidistributed;  then 

(3.5)  X(t*)  -  ?C(t^)  1  ^(t*)  -  ^(t^)  +  c(t*-t^). 
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Proof:   Suppose  that  for  a  particular  n,   a 

n 

satisfies 

a  <  ^ 
n    r 

where   a   is  a  lower  bound  for   X,  ;   then  the   (n+l)^^ 

segment  of  X^      starts  at  distance   r   to  the  right  of  the 

th 
■beginning  of  the  n   segment.   Let  n-,   and  n^  be 

1  d. 

integers  such  that   t^  =  n-^^s^^  and   t*  =  n  s    denote  by 
S  =  S(n-j_,n2,a)   the  number  of  a^,  "-i  1  ^  1  ^2  ^'^^   which 

(3.6)  a   <  ^ 

n    r 


It  follows  that  X^     contains  at  least   S   jumps  to 
ther  right,  at  most  n^  -  n-^  -  S   jumps  to  the  left,  and 
so 

(3.7)  rj^(2S-(n2-n^))  <  X^i^*)    -  ^^{tj  . 
Equidistribution  means  that 

(3.8)  — §—  =  [l+o(l)](l  +  M)/2 
ng-n^  ^ 

which  implies  that 

ri^[2S-(n2-n-j^)]  =  [a+o(  l)  ]  s^(  n^-n^)  . 

Substituting  this  into  (3-7)  we  get  that 


4^5 


a(t2-t^)  +  0(1)  <  7(^{t^)    -  Xi,it^)     . 

Letting  k   tend  to  infinity  we  obtain  that 
a(t2-t^)  <  X(t2)  -  Xi^^), 

which  implies  the  left  side  of  inequality  (3-3);  the  right 
side  follows  similarly,  and  so  does  part  b)  of  lemma  3.2. 

It  follows  from  the  weak  law  of  large  numbers  that 
almost  all  sequences  ^a  ?  are  equidistributed .   In  an 
actual  calculation  we  need  not  rely  on  this  property  of 
random  sequences  but  would  preferably  employ  a  particularly 
evenly  distributed  sequence,  such  as  the  fractional  part 
of   n  9,   where   6   is  a  quadratic  irrationality. 

In  what  follows  we  assume  that  the  sequence  Ca  ? 
is  equidistributed. 

Lemma  3  «3 :   Let   X  -   li'^  Xi,   t)e  a  characteristic 
curve.   Then  we  can  select  a  subsequence  for  which 

(3.9)  lim  Xj^(t)  =  X(t) 

for  all  but  a  countable  set  of   t. 

Proof;   According  to  theorem  2.1,  part  (D),  the 
total  variation  of  Xv,   over  a  given   t   interval  is 
uniformly  bounded  for  all   k.   Also,  the  functions 
X,  (t)   are  uniformly  bounded.   It  follows  from  these 
facts  by  Kelly's  theorem  that  for  a  suitable  subsequence 
the  limits 


iii^ 


lim   ;^,  (t)  =   s(t) 


k 
and 

lim J  WX^{^)  I  -  s(t) 
o 

exist  for  every  t.   We  claim  that  (3-9)  holds  at  every 

point  of  continuity   t   of  's(t).    To  see  this  choose, 

given  any  positive   5,   t„   and   t*   close   enough  to 


t    so  that 
o 


and 


t„  <  t   <  t* 
*    o  — 


r 

t^ 


dXj  <  6 


Clearly  assumption  {^.^).     of  lemma  3.2  is  fulfilled  with 
a  -  s(t^)  -  5,   b  =  s(t^)  +  6: 

s{t^)  -  5  <  X^   <    s(t^  +  6). 

So  we  conclude  by  that  lemma  that 

s(tQ)  -  6  <     X(t)  <  s(t^)  +  &   . 

From  these  bracket ings  of  ^,       and  X    ^^^  ^^e  arbitrariness 
of   6   relation  (3-9)  follows. 

Having  shown  that  the  slope  of  every  characteristic 
X    is  the  limit  of  the  slopes  of  the  approximate  character- 
istics  Xv   of  which  X    is  the  limit,  we  turn  to  showing 
that  likewise  the  value  of  the  solution  u  on  either  side 
of  X  is  "ths  limit  of  the  values  of  the  approximate 
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solutions   u,   on  either  side  of   Xi,-   The  key  fact  here 

is  the  one-sided  equicontinuity  of  the  approximate  solutions: 

Lemma  3 -^ :  We  can  select  a  subsequence  such  that 
the  following  statement  is  true  for  all  but  a  countable 
set  of  values   t: 

Given  any  positive   e   there  is  a   5  =  5(e,t)   such 
that  for  all   0  <  y  <  5 

(3.10)  |u^(;^^(t)  +  y,t)  -  Uj^(;^j^(t)  +0,t)|  <  e. 

Before  presenting  the  proof  of  lemma  3-^  we  draw 
some  conclusions  from  it: 

Corollary:   For  all  but  the  aforementioned  countable 
set  of  t 

(3.11)  lim  Uj^(Xk(t)  +  0,t)  =  u(;t(t)  +  0,t). 

Proof  of  corollary:   According  to  part  (c)  of 
theorem  2.1  the  total  variation   of  u,  (/C,  (t)  +  0,t)   is 
uniformly  bounded.   So  by  Helly's  theorem  we  can  select 
a  subsequence  for  which  the  limit  on  the  left  of  (3.11) 
exists;  denote  the  value  of  this  limit  as  U— (t) .   We 
want  to  show  that  ui(t)  =  u(x(t)  +  0,t). 

By  the  triangle  inequality 

(3.12)  |u±(t)  -  u(7((t)  +  y,t)l  <  |ui(t)  -  Uj^(Xi^(t)  +  0,t)l  + 

+  l^k(Xk(^)  ±  0'"^)  -  Uk(?f(*)  ±  y^t)  I  + 
+  |Uk^X(t)  ±  y>t)   -  u(X(t)  +  y,t) I . 
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The  first  term  on  the  right  is  by  definition  of  U   less 
than   e   for   k   large  enough;  the  second  term  is  by  (3. 10), 
less   e  provided  that   y  <  5/2   and   k   is  so  large  that 
■  IX-^i'^)    -  X("t)  I  <  y.   Finally,  according  to  the  convergence 
theorems  of  [    ]   quoted  at  the  outset  of  this  section 
the  integral  of  the  last  term  with  respect  to  y  over 
6/4  _<  y  ^  6/2   tends  to  zero  as   k  — ^  00  ;   this  implies 
that  there  is  a  sequence   y^^,   in   (5/4,  6/2)   such  that 

|u^(X(t)  +  y^,t)  -  u(X{t)  +  y^,t)| 

tends  to  zero.   So  for   k   large  enough  the  above  quantity 
is  less  than   e;   applying  this  to  the  last  term  on  the 
right  in  (3-12)  with  j  =  y^     we  obtain 

|U-(t)  -  u(}|^(t)  +  yj^,t)|  <  3e. 

Since   e   Is  arbitrary  and  y,   can  be  taken  as  close  to 
zero  as  desirable,  by  choosing   5   small  enough,  it  follows 
that  U-(t)  =  \i{-^{t)    +   0,t).   This  completes  the  proof  of  (3.11) 

Denote  as  usual  u(X(t)  +  0,t)  -  u(X(t)  -  0,t) 
by   [u]. 

Theorem  3.1:   Along  each  characteristic  X  ,   the 
Rankine-Hugoniot  relations 

(3.13)       X(t)  =Iil  =  isl 

[u]    [v] 

are  satisfied  at  all  points  where   [u],  [v]  /  0,   and 
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the  characteristic  relation 

(3.1^)  x(t)  =  p(ua(t),t)) 

at  all  points  where   [u]  =0,   again  with  the  possible 
exception  of  a  denumerable  set  of   t. 

Proof:   By  construction, the  slopes  of  the  approximate 
characteristics  satisfy   the  Rankine-Hugoniot ,  respectively 
characteristic  relations.   Letting  k  tend  to   co   and 
using  the  limit  relations  O .9)    and  (3.11)  we  obtain 
(3.13)  and  (3-1^)  as  limits  of  (3-13)^^   and  (3-1^)^ 
respectively. 

We  turn  now  to  the  proof  of  lemma  '3.^.      We  recall 
from  section  2,  equation  (2.9)  and  (2.10),  the  quantities 

(3.15)         C^ight^^)'  ^left(^)'  Q^^) 

denoting  the  amount  of  right  and  left  wave  cancelled  in 
the  diamond  A,   and  the  amount  of  wave  interaction  in  A; 
these  quantities  depend  of  course  on  the  approximate 
solution  u^.   We  denote  by   dC^^S  ^^[  ^-^^   dQ^^ 

the  measures  which  assign  the  quantities  (3.15)  "to  the 
center  of  A. 

We  further  recall  from  section  2,  equation  (2.22) 
and  (2.l8),  the  quantities 

(3.15)^        xi(X^A)   and   S(A) 

denoting   the   ainount   of   left  wave   crossing     ^    in     A,      and 
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the  amount  of  right  shock  entering   /^  in  A.   We  denote 
by  dX^  y/Yland   dS   /^the  measure  which  assigns  the 
quantities  (3-15)    to  the  center  of  each  A   intersected 
t)y  /^  ;  the  superscript   k   calls  attention  to  the  dependence 
of  these  measures  on  u,  . 

We  select  a  subsequence  so  that   ^ v  — ^  X  ^^^    ^° 
that  the  limits 

(3.16)        dC^^^  -^  dC^,  dc]^^  -^  dC^jdQ^^^  -^   dQ 


(3.17)      dx±^(Xk)  -^  dxi(X),  ds^^^Xk)  -^  ^s(X) 


exist  in  the  w*  topology  for  measures. 

Let   t   be  any  value  of   t   for  which  the  point 
X(t  ),  t   has  measure  zero  with  respect  to 

dC^  +  dC_g  +  dQ 

and  such  that   t   has  measure  zero  with  respect  to 

o 


dX+(X)  +  |dX"(X)l  +  dS(X) 


The  excluded   t   form  at  most  a  countable  set. 

o 

Let  7  be  any  sufficiently  small  positive  number; 
we  shall  denote  by  0(7  )   any  quantity  less  than  K7  , 
K   independent  of  7   and   k.   For   t   as  above  we  can 
find  a  neighborhood  Y{j)      of    ^(t  ),  t    such  that  for 
all  k  large  enough 
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(3.18)  The  a-mount  of  wave  In  u,   cancelled  in  V   is 
less  than  j'^ 

(3.19)  The  amount  of  left  wave  crossing  7^,       in  V   is 
less  than  y^ 

(3.20)  The  amount  of  shock  entering   X 1^  }^   V  is  less 
than  7 

(3.21)  The  amount  of  interaction  in  V  is  less  than  7 
We  shall  apply  the  conservation  laws  and  derive 

some  further  estimates  on  waves  in  V.   The  final  upshot 
will  be  that  the  oscillation  of  u,   at   t   on  either 

K  O 

side  of  X,       along  an  interval  whose  length  depends  only 
on  7   is  small  if  7   is  small.   This  is  inequality  (3-10) 
of  lemma,  3  -^  • 

Consider  the  domain  vv  ,   shown  on  figure  3-l> 
hounded  "by   Xi  ,      ^^   interval  cr    located  at   t  <  t 

K  O 

and  some  time-like  line   t,   all  located  in  Y: 


t  A 


t   + 

o 


Fig.  3.1 

We  apply  the  law  of  conservation  of  left  waves  to  -A.  . 
Since  (3.18),  (3. 19)  and  (3.21)  estimate  the  amount  of 
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cancellation,  the  amount  of  left  wave  crossing  y( ,  , 
and  the  amount  of  interaction,  respectively,  we  obtain, 
see  (2.25),  an  estimate  on  the  amount  of  left  wave 
crossing   o*  : 


(3.22)  l^left^^^l  -  ^^^  • 

Using  once  more  the  law  of  conservation  of  left  waves 
for  a  domain  bounded  by  any  right  characteristic   0, 
contained  in  a  subdomain  of  V   (to  be  again  denoted  by 
V) ,  a  spacelike  line  cr     as  above  and  a  timelike  line, 
and  using  (3.I8)  and  (3.21)  we  deduce  as  before  that 

(3.23)  l^left^^^l  -  0^>'^^- 

(3.22)  is  an  estimate  for  the  total  amount  of  left  waves 
crossing  c  ;   we  wish  to  estimate  also  the  total  amount 
of  right  waves  crossing  cr.      We  start  with  proving  the 
following 

Assertion;  Let   0,   be  any  approximate  characteristic 
located  in  V  which  contains  at  time   t   a  shock  of 

strength   >  y   '     .   Then  the  strength  of  0,   at  any 

3/2    3 
later  time   t'  >  t   is  at  least  y^'       -   j    ,      at  most 

Proof;   According  to  (3.I8)  the  amount  of  shock 
cancelled  in  V  is  at  most  7-^;   by  (2.27)  "the  strength 
of  0.,   at  times  before  it  has  attained  the  strength 
y    '        is  at  most  y   '      +  7  ,   and  afterwards  at  least 
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7  -  7  •  To  estimate  the  subsequent  strength  of  0, 
from  above  we  denote  by  |_  P J  the  shocks  which  enter 
0,  .   Since  the  strength  of  0,   is  at  least    p/    > 

the  entry  of  a  shock  of  strength  P   causes  an  interaction 

1  3  /2 
in  the  amount  at  least  i<y        P.   Since  by  (3.21)  the  total 

Interaction   in  V  is  at  most  y    ,      by  (2.28) 

(3.24)  y-  p  <  27^/2. 


According  to  the  discussion  in  section  2,  the 
paragraph  preceeding  relation  (2.23),  the  variation  of 
u,   along  either  side  of  a  characteristic  0,   is 
bounded  by  the  sum  of  the  total  amount  of  shock  entering 
0,  ,   the  amount  of  right  wave  cancelled,  and  the  amount 
of  left  wave  crossing  0,  .   These  quantities  are  estimated 
by  (3.24),  (3.18)  and  (3-23).   According  to  (2.23)  the 
total  variation  of  u,   on  either  side  of  0,   furnishes 
an  estimate  on  the  total  variation  of  the  speed  ^ 
of  0,  .   Thus 

( 3 . 25 )  The  total  variation  of  the  speed  of  any  right 
characteristic  0,   in  V  which  starts  with  a  shock  of 
strength  >  j'' '        is  at  most  0{y   '     )  . 

To  proceed  further  we  have  to  distinguish  two  cases, 
corresponding  roughly  to  the  limit  characteristic  "yC 
being  a  shock  of  strength  _>  7   or  not. 

A)   There  is  a  right  shock  a   contained  in  u, 
of  strength   >  7   for  all  k  large  enough,  and   a 
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tends  to   y(t  ) ,  t  . 

B)   V  can  be  chosen  so  small  that   for  infinitely 
many  k   all  right  shocks  of  u,   in  V  have  strength 
less  than  7. 

We  choose  R  =  R(7)   so  small  that  a  square  centered 

at   y(t  ),  t   with  edge   12R   is  contained  in  V. 
'*•  ^  o     o         ° 

We  shall  refer  to   12R  as  the  inner  diameter  of  V. 

We  turn  now  to  case  A) ;  we  denote  by  f,       that 

portion  of  the  characteristic  through  a        which  is 

contained  in  V.   Assuming  that  all  characteristic  speeds 

are   <  1   (easily  accomplished  by  rescaling)  we  conclude 

that  for   t   <  t  <  t   +  6r,  ip,       lies  in  V. 
o  —   —  o     '    ^K 

By  assumption  A)   for  k   large  enough  the  initial 
point  of  ij/-.       lies  close  to   X(t  ),  t  .   We  shall 
eventually  show  that  f,       in  fact  lies  on  X,  . 
As  a  first  step  we  show  that  all  characteristics  point 
toward  tp-.;      more  precisely: 

Assertion;   There  exists  a  positive  constant   c 
such  that 

(3.26}  p(Uj^(r,t))  +  C7  <  i^{t) 

for  all   t   satisfying   |t-t  ]  <  R  and  all   r   such  that 

(3.27)  ^^(t)  <  r  <  ^^(t)  +  C7R  . 

Similar  estimates  hold  to  the  left  of  tp,. 
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Proof  of  assertion:   Since  the  strength  of  Tp, 
initially  is  _^  7?   it  follows  "by  our  previous  reasoning 
that  its  strength  is  always  >  y    -   y    .       It  follows  from 
the  genuine  nonlinearity  of  the  system  under  study,  see 
section  1,  equation  (1.4),  that  the  characteristic 
speeds  of  the  state   u   ^    on  the  right  side  of 
f,       satisfies 

(3.28)  pCu^^S'^*)  +  c^7  <  i^^    . 

This  shows  that  (3.26)  holds  for   r   adjacent  to  i^... 
We  shall  prove  that  (3.26)  holds  with  c  =  c  /2   for  all 
r   satisfying  (3.27).   For  suppose  on  the  contrary  that 

(3.28)  is  violated  at   t  =  t,  ,  l^^-t  |  <  R,   at  some 
point,  say,  to  the  right  of  f,       satisfying  (3.27).   We 
abbreviate   ^,  (t,  )   by  p,  ,   and  denote  by   r,   the  first 
point  to  the  right  of  p,   where  (3.26)  is  violated  at 
time   t,  .   We  claim: 

There  is  a  right  characteristic  through   (r,  ,t,  ) 
whose  speed   s   satisfies 

(3.29)  s  =  i^^i\)    -    C7  +  0(7^). 

Proof:   According  to  (3.28),  (3.26)  holds  at  the  point 
adjacent  to  tp-.;      therefore  (3.26)  is  satisfied  at  the 
constant  state  on  the  left  of  the  waves  issuing  from 
r,  ,   while  it  is  violated  at  the  state  on  the  right. 
Those  two  states  are  separated  by  a  left  wave  and  a  right 
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wave;  according  to  (3.22)  the  strength  of  the  left  wave  is 
Oiy    ) ,      and  therefore  across  it   p   changes  by   0(7  ). 
Since   p   decreases  across  a  right  shock  wave  and  increases 
across  a  right  rarefaction  wave,  it  follows  that  the  right 
wave  through  r,   is  either  a,  shock  wave  of  strength 
0(7  )   or  a  rarefaction  wave;  in  either  case  (3.21)  can 
be  satisfied. 

Denote  by  0,       the  right  characteristic  through 
(r,  ,t,  )   with  initial  slope   s,   restricted  to  V: 


t  A 


\-L 


^k 


0^. 


Fig.  3.2 


Then  (3.29)  can  be  rewritten  as 


(3.30) 


^^(t,  )  =  f,  (t,  )  -  C7  +  0(7^) 


'k^  "k- 


^k^"k 


Denote  by  <y,       the  interval  p,  <  x  <  r,,    t  =  t,  ; 


k 


k  -  ^  -  k^ 


^k 


We  claim: 


(3-31)   All  right  shocks  crossing  the  segment  cr,       have 
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strength   <  y' '         for  k   large  enough. 

Proof:   Suppose  on  the  contrary  that  some  right  shock 
p,   has  strength   >  y  fo^  arbitrarily  large  k.   Denote 

by   K,   the  right  characteristic  issuing  through  p,  , 
restricted  to  V.   Since  p,   lies  between  p,   and 
r,  ,   its  speed  satisfies  (3 .26): 

« 

According  to  (3.25)  the  speed  of   rc   and  of  ^,   varies 

3/2 
at  most  by  0{y  '     ),  and  so 

(3.32)  Kj^(t)  <  i/^(t)  -  C7  +  0(7^/2) 

for  all   t. 

By  contrary  assumption  t,  -  t   <  R;   R  was  so  chosen 

that  the  square  with  sidelength   12R   centered  at 

'^(t  ),  t    is  contained  in  V.  Since  we  have  assumed 
'  o  '   o 

that  all  characteristic  speeds  are   <  1   it  follows  easily 

that  the  portion  t,  <  t  <  t,  +  2R  of  both  ^,   and   K 

are  contained  in  V. 

According  to  inequality  (3.32)  the  speed  of   (6   is 

less  than  the  speed  of   ^,   by  Q.y   +  0(7  )  .   Therefore 

according  to  part  b)  of  lemma  3 -2^ for  k   large  enough 

tp-.       and  !<■,       approach  each  other  at  the  rate  of   C7  +  0(7  )  . 

On  the  other  hand,  the  distance  of  tli,       and   /C,   at   t  =  t, 

'  ^k        k  k 

is  according  to  (3-27)  less  than  C7R;  therefore  the  distance 
of   ^,   to  ic       at   t  =  t,  +  2R   is  less  than 
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cyR  -  ic-y+0{-y^)2R     which  is  negative  for   7   small  enough. 
Thus  Tp-.       and   ic,   meet  at  some  time  between  t,   and 

t,  +  2R;   since  the  first  shock  has  strength  y,      the 

3/2 
second  y        ,      in  the  diamond  where  they  meet  there  is 

5/2 
interaction  in  the  amount  of  y^'    ,      contrary  to  (3.21). 

From  the  estimate  (3 -31)  we  deduce 

(3.33)   The  total  variation  of   z ( x , t )   along   cr^   Is^ 
0(7^). 


k 


Proof:   The  Riemann  invariant   z   changes   across 
left  waves  and  right  shock  waves;  according  to  (3.22) 
the  total  amount  of  the  former  is  bounded  by   0(7  ). 
Across  a  right  shock  of  strength  P   the  change  in   z 
is   0(P  );   therefore  the  total  variation  of   z   across 
right  shocks  is  bounded  by 

MaxP^ 

According  to  (3.31)   MaxP   over  c,       is  y^'     i      as  shown 
in  Section  2,  ^^  p.   is  bounded  a  priori  by  the  total 
variation  of  the  initial  data.   Thus   Maxp  YZ   ^  •  1. 
[_Y       )   const;   this  completes  the  proof  of  (3.33)- 

Next  we  show  that  not  only  the  shocks  crossing  the 
interval   cr:   but  all  those  contained  in  V  and  lying 

between   (//,   and  0,   have  for   t  >  t,   strength  less 

3/2 
than  Q){y        ).   To  show  this  we  denote  by  P    the  strength 
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th 
of  the  strongest  right  shock  wave  at  the  n   time  step 

which  lies  between  ip,       and  0,.      According  to  the  "basic 

theorem  governing  the  interaction  of  waves,  see  equations 

(2.7),  (2.8),  two  shock  waves  of  strength  at  most   P 

can  produce  a  shock  wave  of  strength  at  most   2P   plus 

an  error  term.   Thus 


(3.3^)  p  ,,  <  2p   +  Q 

\^  >  /  "^n+l  —  ^n   ^n 

where   Q   is  the  amount  of  interaction  in  the  diamond 
from  which  3   ,   originates. 

It  is  easy  to  see  that  in  the  production  of  a  shock 
of  strength  p  ^  out  of  two  shocks  of  strength  at  most 
P  ,   the  amount   Q   of  interaction  satisfies  the  inequality 


"^n 


(3.35)  (^n+l-V^n^^VQn^  - 


Adding  ^^   -   4q^  to  (3-3^)  we  obtain 


^n+1   ^n    ^n  -  -^^^n  ^n'' 


using  this  to  estimate  the  second  factor  on  the  left  in 

(3.35)  we  get 

(3.36)  (Pn+i-V^n^Pn+l+V^Qn)  -  5^   ' 


Denote  the  maximum  of  P   by  B;   let  N  be  the  smallest 
index  for  which  p   =  B.   We  sum  (3-36)  over  n,  ^  n  _<  N; 
using  the  identity 
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Vl-^n)(Pn+l+^n)  =  ^n+l'^n' 


the  estimate  (3.21) 


Qn->'^^ 


the  definition  of  B: 


M  B,  Pj,  =  B 


and  estimate  (3-31) 


we  obtain  that 


B 


2    „3 


7^  -  6B7^  <  0(7^) 


which  clearly  implies  that 

(3.37)  B  <  0(7^/^), 

as  asserted. 

Next  we  claim  that 

(3.38)  ^j^(t)  <  ^^(t)  -  C7  +  0(7^^^) 

for  all   t  >  t,  ,   k  large  enough. 

Proof;   According  to  (3.30),  (3.38)  is  true  at   t  =  t,  . 
According  to  section  2,    the  changes  in  the  direction  of 
a  right  characteristic   0,   are  due  to  crossings  by  left 
waves,  cancellations  and  entering  shock  waves.   The 
contributions  of  the  first  two  are  according  to  (3.23) 
and  (3.18)   0(7  );   to  estimate  the  effect  of  shocks 
entering  0,       we  distinguish  shocks  entering  from  the 
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right  and  those  entering  from  the  left.   Those  entering 
0,   from  the  right  decrease  the  speed  of   0,   and  therefore 
favor  (3.38);  so  to  complete  the  proof  of  (3-38)  we   have  to 
show  that  the  total  amount  of  wave  which  enters   0,   from 


3/2 
the  left  does  not  exceed   0(7  '^  ) 


Denote  by  P    the  strength  of  the  shock  wave  which 

th 
enters   0,   from  the  left  at  the  n   time  step;  since  the 

strength  of  0,   at  that  time  is  at  least  ^^  P   -  C  . 

m<n 
where   C    is  the  amount  of  wave  cancelled  in  0,   up  to 

time   n,   the  quantity 


^n^'^   ^m  n' 
m<n 

th 
is  a  lower  bound  for  the  interaction  at   0,   at  the  n 

k 

time  step.   Summing  with  respect  to   n  and  using  the 
bound  (3.21)  for  the  total  amount  of  interaction  we  get 

;^pp     -"5~Cp  <7^. 
m<n 

This  can  be  rewritten  as 

1 


^^^—Txi'  2  ^—   ^n   ^—     XX  n  —  ' 


According  to(3.i8),  C  <   -p    and  according  to  (3  •37)  shocks 

n  —  ' 

to  the  left  of  0,   satisfy  p   <  0(7^/^).  Using  these 

K  n  — 

estimates  in  the  above  inequality,  and  abbreviating 

P    as  'J~     we  obtain 
n      ^— 
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1  ^-2    ^,    3/2  V  ^-    3 
from  this  we  deduce  that 

as  asserted. 

Inequality  (3-38)  relates  the  speeds  of  f.       and 
0j^;  it  follows  then  from  part  t> )  of  lemma  3.2  that  ip      and  0 
approach  each  other  at  the  rate  cj  +  Qi{'y    '    )  .   Since 
according  to  (3-27)  at  time   t    they  are  less  than 
C7R  apart,  it  follows  as  before  that  at  some  time 
^  "  ^cross   between   t^^  and   t^^  +  2R,  -^^  and   0^   run 
into  each  other.   We  apply  now  the  law  of  conservation 
of  right  rarefaction  waves  to  the  region  vA.  bounded  by 
a     and  by  the  segments  of  -^^     and   0,   between   t,   and 
"^n^n.ci^'        No  right  rarefaction  wave  leaves  -A.  across 
^k  ^^     ^k*   therefore  we  conclude  that  the  amount  of 
right  rarefaction  wave  crossing   c,  is  equal  to 
C^^^(-A.)  +  Q(^^(y\u)   which  by  (3.18)  and  (3-21)  is   <  27^. 

Since  the  only  kind  of  wave  across  which  the  Riemann 
invariant   w   increases  is  a  right  rarefaction  wave,  it 
follows  from  what  we  have  Just  shown  that  the  increase 
of   w   across   o-^  is  at  most   0(7^).   According  to  (3.22) 
the  change  in   z   across  <y     is   0(7^);   since  according 
to   (1.4)  p   is  an  increasing  function  of  w   it  follows 
then  that  the  increase  of   p   across   cr  is   0(7^)  .   On 
the  other  hand  combining  (3.28)  and  (3.29)  we  get 
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p^^right^  <  p(u(r,  ,t,  )  -  C7  +  0(7^) 


k'  k' 

which  indicates  an  increase  in  p  across  cr  by  cy  +  0(7  ). 
This  is  a  contradiction,  into  which  we  got  by  denying  (3.26), 
thus  proving  the  truth  of  (3.26). 

By  hypothesis  A)  the  beginning  of  -^^   lies  close  to 
X(t  ),  t   and  so  inequality  (3.26)  is  applicable  to  the 
states  on  either  side  of  Xv *   ^^  conclude  from  that 
inequality  and  lemma  3.2b)  that  ^,   and   ^',   approach 
each  other  at  the  rate  of  07   and, since  initially  they 
are  not  far  apart  they  must  meet   in   V  for   k   large 
enough.   When  they  meet  a  shock  of  strength  7   would 
enter  X.,      contradicting  assumption  (3.20) .   The  only 
way  out  of  this  contradiction  is  that  il/,       itself  lies 

on  Xj^. 

Applying  (3-26)  to  Xv  we  conclude  that  all  character- 
istics within  a  distance  of   C7R  of  .^,  ,   and  for 
|t-t  I  <  R,   approach  Xv^   at  the  rate  of   07.   Thus  all 

characteristics   0,   which  originate  at   t  =  t    inside 

'^k  *="  o 

the  interval 

(X(tQ)-C7R,X(t^)+C7R) 

must  for  k  large  enough  run  into  Xi,-   Applying  the  law 

of  conservation  of  right  shock  waves  to  the  region  bounded 

by  X,  ,  0,   and  an  interval   cr  on   t  =  t    we  conclude 
^    k'  '^k  o 

that  all  but   27   of  the  amount  of  a  right  shock  crossing 
(J     must  enter  either  Xi,   o^  ^1,  •   By  (3-20)  no  more  than 
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7    can  enter  ^,  ;   on  the  other  hand  the  strength  of  (^ 

2 
cannot  exceed   7    since  then  there  would  be  more  than 

7   amount  of  interaction  with  /j^  when   0j^   runs  into 

%,,      contrary  to  (3.21).   This  shows  that  the  amount  of 

2      3 
right  shock  crossing   C  is  at  most   7   +  0(7  ). 

Finally  we  deduce  from  the  conservation  of  right 

rarefaction  waves  that  the  total  amount  of  right  rarefaction 

wave  crossing  Cf     is  0(7  ).   Since  the  total  variation  of  w 

and   z   along   c  can  be  estimated  in  terms  of  the  total 

amount  of  wave   crossing,  we  deduce  that  the  total  variation 

of  w   and   z   on  either  side  of  /t,   at   t  =  t    is  leas 

k  o 

2       3 
than  7   +  0(7  );  but  then  so  is  the  total  oscillation. 

Since   7  may  be  taken  arbitrarily  small,  this  completes 

the  proof  of  lemma  3-^  in  case  A. 

We  turn  now  to  case  B;  here  by  assumption  all 
right  shocks  are  weaker  than  7.   It  follows  from  this 
and  from  estimate  (3.21)  for  the  interactions  that  the 
amount  of  shock  entering  any  right  characteristic  is  at 
most   0(7).   Using  the  arguments  leading  to  the  proof  of 
estimate  (3.25)  we  conclude: 

The  total  variation  of  the  speed   ^,   of  any 
right  characteristic  in  V   is  less  than 

(3.39)  K7 

K   independent  of  k  and   7. 
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We  use  this  to  prove  that  all  characteristics  point  away 
from  each  other.   More  precisely: 

Assertion:   Choose  as  before   R  =  RCy)   so  small 
that  a  square  of  side  length  12R   centered  at  "/(t  ),  t 


o 


is  contained  in  V.   Let   (a,t-,)   and   (b,t-,)   denote  two 

points  in  V,   with  t-,  <  t   and 
^  '         1  —  o 


(3. 40)  ■     X(^i)    <  a  <  b  <  X(t^)  +  6KrR 

let  0^  and   0^    denote  approximate  right  character- 

istics through  these  points;  then 

(3.41)  ^l^Ut^)    <  9^l^Ut^)   +   3K7. 

A  similar  inequality  holds  to  the  left  of   X • 
Proof:   Suppose  (3 .41)  is  violated;  then,  since 

t>y  (3-39)  characteristic  speeds  oscillate  only  Ky , 

we  find  that  for  all   t  >^  t  , 

(3.42)  '0l^Ut)    >   9il^ht)    +  K7  . 

According  to  lemma  3-2,  part  b),  this  implies  that  for 
k  large  enough  0^    and   0^  '      approach  each  other  at 
the  rate  Ky .   Since  at  time   t-,   they  are  less  than 
6K7R  apart,  it  follows  that  they  must  intersect  at  time 
less  than  t-,  +  6r  .   Assuming  that  all  characteristic 
speeds  are   <  1,   the  intersection  lies  in  V.   We  apply 
now  the  conservation  law  for  right  waves  to  the  region 
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/»,   bounded  by  0l_      ,    0l_    '      and  the  interval   (a,b) 


at   t  =  t 


-,  .   We  conclude  that  the  amount  right  rarefaction 


wave  crossing   (a,b)   is  bounded  by  C^  ^^k^  "*"  ^  ^k^ 
which  according  to  (3 .18)  and  (3.21)  is   <  2y^ .      The 
amount  of  right  shock  wave  crossing   (a,b)   is  bounded 
by  2-y        plus  the  amount  which  enters  the  bounding 
characteristics  0^  or  01      •      Since  the  strength  of 

these  shocks  cannot  exceed  7,   and  the  amount  of 

■5 
cancellation  on  them  is  at  most  7'^,   the  amount  entering 

them  cannot  exceed  2-/   +  7  ,   and  so  the  amount  of  right 

shock  crossing   (a,b)   is  at  most  2-y    +  37    -      Combining 

this  with  (3.22)  according  to  which  the  total  amount  of 

left  wave  crossing   (a,b)   is   0(7  ),   we  conclude  that 

the  total  variation  of  u,   in   (a,b)   is   0(7);   in 

particular  then  the  speed  of  the  characteristics  at   a 

and  b   differ  by  less  than  Kj ,      contradicting  (3.42) 

and  proving  our  assertion. 

We  turn  now  to  deriving  an  inequality  in  the 

direction  opposite  to  {3.^1).      For  this  purpose  we  look 

at  the  characteristics  01^        issuing  from  the  points 


(3.43)  ^(^o-^)  +  y,  t^  -  R 
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o 


^o  -  R 


Fig.  3-3 

There  may  be  actually  a  whole  fan  of  characteristics 
issuing  from  point  (3  •^3). 


■(y) 


'(y) 


We  denote  by  J^    the  tangent  line  to  0^^-^ 
at  the  point  (3-'^3)- 

Using  (3.^1)  with  a  =  X^(t^-R),  b  =  :^^(t^-R)  +  y 
we  conclude  that 


By  (3-39)   Xi,   oscillates  less  than  Ky ;      using  lemma  3-2 
part  b)  we  conclude  from  (3-^^)  that  for   K  large  enough 


^k(to)  -^k(to-R)  -KrRlli^^t  )  -fl^Ut-n)   +3K7R 


■k   ^  o'    -J-k   ^  o 


Since  ^il"^  (t  -R)  =  Xy-i'^    -R)  +  y»   ^^  conclude  that 


(3.45)  ^k(to)  -  4K7R  +  y  <  $jj.^^t^)   . 

On  the  other  hand,  since  0i'^        is  X,  ,   we  conclude, 
again  using  (3.39)j  that 

(3.46)  Ik°^(^o^  -  ^k^^o^  +  K^^   • 


¥e  claim  that  there  is  a  value  y,  , 


(3.4?)  0  I  yk   ^  ^^^^ 

such  that 

(3.48)  Ik^^^t^)  =Xy.{t^)    +   2K7R  +  0(7^R). 

Proof:   As   y   increases  $il   ( t  )   jumps  to  the 
right  only  across  left  shocks,  whose  magnitude  according 
to  (3.22)  is   0(7^);   and  by  (3.^5)  and  (3-46)  the  value 
of  "^  Z        a-t  the  left  and  right  endpoint   of  the  interval 
(3.47)  is  less,  respectively  greater  than  the  right  side 
of  (3.48). 

Denote   0,      as   ?//,  ;   we  show  that  Tp     is  distinct 
from  "X    but  nearly  parallel  to  it.   Applying  once  more 
(3.39)  we  conclude  that  for   k   large  enough  ^i,(t  ) 
differs  from  Y^(t  )    at  most  by  K7R;  so  from  (3.48) 
we  conclude  that,  mod  0(7  R) 

(3.49)  ^(^o^  +  ^^^  "  ^k^^o^  "\^^o^    +  ^^^^ 
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We  use  the  right  portion  of  this  inequality  to 
estimate  the  average  speed  of  ip,;      using  also  the  fact 
that  Tp,       lies  to  the  right  of  Xu.     we  get 


•^k^  o'    "^k^  o  "^'  -'"-k^  o'    ^'^'''        ''K'    o  ''" 

^k(^o^  -  ^k(^o-^^   ^k^^o)  -  ^k(^o-^^  ^  ,^, 
<  +  3K7. 


R  R 

Using  (3.39)  and  lemma  3.2  we  conclude  that  for  k   large 
enough  the  average  speed  of  a  characteristic  differs  from 
its  speed  at  any  point  at  most  by   K7;   so  we  conclude 
that 

(3.50)  ^k^^o^  -  >^k(^o^  +  5Kr 

Now  suppose  that   a  lies  in  the  range 

(3.51)  Xk(to)  <   a  <  ;j:i,(t^)  +  K7R; 

according  to  the  first  inequality  in  (3. ■'+9),   ^^^   lies 

to  the  right  of  0^^^-   Applying  (3.^1)  to  0^^^   and 

to   "1^11=  ■^v  (with  t,  =  t^)   we  conclude,  using  (3.50), 


that 


^k^^(^o^  -  ^k^^o^  +  ^^^  -^k^^o^  +  ^^^ 


On  the  other  hand,  since  0l^   lies  to  the  right  of  ^,  , 
we  can  apply  (3.^1)  to  deduce 
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Thus  we  have  upper  and  lower  "bounds  on  0!_    for  all  a   in 
the  range  (3-51)  • 

Since  the  strength  of  shocks  is   <  j,    0,   differs  from 
p(u,  )   at  most  by   0(7)   and  so  we  can  conclude  that  the 
oscillation  of   p ( u  )   in  the  interval  ( 3 • 5 1 )  i^  0(7). 

According  to  (3-22)  the  amount  of  left  wave  crossing  (3 -51) 
is   0(7^);   since  right  shocks  are  weaker  than  7  we  conclude 
that  the  total  variation  -  and  a  fortiori  the  total  oscillation  - 


o 


f  the  left  Riemann  invariant   z   over  (3 -Si)  is   0(7^). 


By  assiimption  the  system  under  study  is  genuinely  nonlinear, 
i.e.   0   7^  0.   Therefore  we  can  estimate  the  total  oscillation  of 
the  other  Riemann  invariant  w   over  (5.51)  in  terms  of  that  of 
z   and   p   and  find  it  to  be  also   0(7).   Since   7  may  be  taken 
arbitrarily  small,  this  completes  the  proof  of  lemma.  J>.k   in 
case   B. 

We  are  now  ready  to  pass  to  the  limit.   As  stated  at  the 
beginning  of  this  section  it  follows  from  the  uniform  bounds 
derived  in  [ 2  ]   on  the  total  variation  along  space-like  curves 

of  the  approximate  solutions   u,   that  a  subsequence  of  them 

loc 
converges  to  some  limit   u   in  the   L-,     topology  in  x, 

uniformly  in  t .   It  follows  from  the  uniform  bounds  obtained 

in  Section  2  on  the  amount  of  right  and  left  wave  cancelled  that 

for  a  further  subsequence  the  measures   dC^    and   dc\    have 

limits   dC  ,  dC'   in  the   w*   topology.   Similarly  the  limit 
3r     Jif 
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dQ  — ^  dQ   exists . 

Let   ?(,  (t)   be  a  sequence  of,  say,  right  approximate 


characteristics,  whose  initial  point   OTv^ '^k  ^ '  ^k   "tends  t 


o 


some  point   x  ,  t  .   According  to  lemma  3«lj  a  subsequence 
tends  to  a  limit   7^{t),   and  according  to  the  corollary  to 
lemma  J).H-   the  limits 

11m  u^{X^{t]+0,t), 

which  we  denote  as   u(X( t )+0, t ) ,   exist  for  all  but  a  denumerable 
set  of  t.   According  to  Theorem  J).l,     Xlt)  =  lim  Xi,!*)   exists 
for  all  unexceptional  values  of   t,  and  is  equal  to  [f ] /[u] =[q] /[v] 
at  points  of  discontinuity,  to   p(u)   at  all  points  of  continuity. 

The  functions   u(X'(  ■fc)+0,  t )   as  defined  above  are  functions 
of  bounded  variation;  since  they  were  defined  in  terms  of  the 
values  of  u,   along   Xv^ 

w(X(t)-0,t)  -  w(X{t)+0,t)  =  |[w]|  =  lim  Strength  X^(t). 

According  to  (2.22)  the  total  amount  of  left  wave  ^(Xu.) 
crossing  Xv  ^^  bounded  uniformly  for  all  k.  Therefore  its 
limit 

X^(X)   and  X"(X) 

exists  as  a  measure  on   7(  • 

Clearly  the  limiting  form  of  Theorem  2.1  holds. 

In  the  estimates  derived  in  Section  4  we  shall  employ  a 
denumerable  set  of  characteristics,  each  the  limit  of  approximate 
characteristics  of  the  same  subsequence  of  the   u  ,   selected  by 
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the  diagonal  process. 

Let   /t-,   and   Xo  ^^  a  pair  of  right  characteristics, 
limits  of   X-|i^   and   "Xov*   Denote  t>y   c,,   approximate  space- 
like curves  constructed  hy  connecting  adjoining   S¥  edges  of 
diamonds,  starting  at  time   t  =  ns   on   ^pi^   ^^^  ending  on  '^■i  i^ 
at  time   t*.   It  follows  from  the  argument  given  in  the  proof  of 
lemma  3*2  that  as   k  tends  to   oo   the  curve   c    tends  to  the 
segment   c,   starting  at  Xp(t),t,   ending  at  %,(t*),t*   travel- 
ling with  speed   -  r/s . 

As  in  Section  2,  we  let 


he  the  amount  of  right  rarefaction,  respectively  shock,  wave 

in  u,   which  crosses   C,,  .   In  the  case  of  a  right  shock  wave 
k  tk  ° 

lying  on  an  endpoint   'Y-i^^nc,,   of  oVv?   we  make  the  convention 
By  a  change  of  independent  variables 


t  =  t,  X  =  X  +  at 

achieve  that 
points 


we  can  achieve  that   A   =  -  P^-   It  follows  then  that  at  all 

o     '^o 


p  =  1a1  +  0(e). 

¥e  can  choose   r/s   to  be   p   +  0(e)   and  satisfy  the  required 
condition 

r/s  >^  Max(  p,  ]  A  j  ) ; 

with  this  choice  the  speed  of  a  space-like  curve  along  which 
t   is  monotonically  increasing  (decreasing)  differs  by   0(e) 
from  the  speed  of  any  right  (left)  characteristic.   This  fact 
will  be  exploited  in  Section  4. 
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that  the  wave  does  not  cross   ^-f-i^* 

Lemma  3 • 3 •   After  passing  to  a  subsequence  the  quantities 

exist  for  each   t   and  have  bounded  total  variation  as  functions 
of  t. 

Proof;   Denote  by  A,     ,  the  strip  bounded  by  ^-i^^  '^Pk^ 

cr.    ,   and   c,  ,  .   Since  the  sides  are  approximate  right  charac- 
teristics,  right  waves  enter  only  through  the  bottom,  rarefaction 
waves  leave  only  through  the  top,  and  shock  waves  leave  through 
the  top  and  possibly  through  the  sides.   ¥e  use  these  facts  in 
conjunction  with  the  conservation  laws  of  Section  2  to  conclude 
that 

H 

We  sum  this  relation  over  any  set  of  disjoint  time  intervals 

(t-,,to)   and  find  that  the  total  variation  of   X— ^  {cr.,  )      is 
^  1^  2  r    ^  tk 

bounded  by 

a  quantity  which  according  to  Theorem  2.1  is  uniformly  bounded 
for  all  k.   Applying  Kelly's  selection  principle  the  conclusion 
of  lemma  3 '5  follows. 

Corollary;   Denote  by  /\     the  strip  contained  between  "X-,  > 
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'v^,  cr^        and  cY,^    ,      where   t-j^  <  t^.      Passing  to  the  limit  in  the 
conservation  law  we  obtain  the  inequality 

(3.53)     1x^(0-^  )|  <  1x^(0-^  )|  +  0{e)Q(A). 

The  family  of  space-like  intervals   c   constructed  above 
can  "be  replaced  by  any  other  monotonic  family,  such  as  the 
horizontal  lines   cr'    obtained  as   lim  cr'  ,   the  union  of 
adjoining  WN   edges  alternating  with  NE   edges.   As  before, 
we  deduce  the  existence  of  the  quantities   X— ( cr,  )   and  the 
boundedness  of  their  variation  as  functions  of  t.   Similarly 
we  deduce  the  analogues  of  the  conservation  laws  (3.53). 

Next  we  wish  to  estimate  the  variation  of  the  amount  of 
left  wave  crossing  cr.       and  so  we  apply  the  law  of  conservation 
of  left  waves  to  A.     .  ,  .   This  time  left  waves  may  enter  also 
through  the  right  and  leave  through  the  left  boundary.   We  obtain 

+  0(e)QX^t2k)  +^  /  1^4^^  H. 

*1 

As  before  we  conclude  that  the  variation  of  Xj{cr,     )    as  function 
of  t   is  uniformly  bounded  for  all   k,   and  so  by  Kelly's 
selection  principle 

X±(o^^)  =  limxj(o'^j^) 

exists  for  all   t   for  a  suitable  subsequence,  and  is  of  bounded 
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variation  as  function  of   t. 

We  need  an  estimate  not  only  for  the  variation  but  also 

for  the  magnitude  of  X— {  o ,  ) ;  for  this  purpose  we  consider  the 

So      ^ 

triangular  region   pV,   formed  by   cr;  ,  Xo\r      ^^"^  '^+\r>      where 


tk 


tk'  '^2k 


■tk^ 


T.,   is  a  time-like  line  starting  at   X-|(t*),  t*    and  ending 
at  some  point  of  Tto : 


Figo  ^.k 


Left  waves  enter   T^vf   only  through   ^ovi   therefore  we  conclude 
from  the  conservation  law  that 


t-t 


k 


where   t-t,   is  the  time  when  t,,  intersects  X^.      Letting  k 
tend  to   00   we  obtain 
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t-t 

o 


We  recall  from  Section  2  that  the  total  arao-unt  of  wave 

(k) 
crossing  a  curve  C     is  denoted  by  Xlc): 

According  to  (2.26),  the  change  of  the  Riemann  invariants   w 
from  one  end  of   C  to  the  other,  denoted  as   [w](o'),   is 

[w^^^K^)  =x:(^^o-)  +X-(^\a)  -f  0(e2)x(^^cr), 


and  similarly  for   z.   Actually  for   z  we  only  need  the  estimate 

|[z(^)](cr)|  <X+(^^(>)  +  lx-(^^a)|  +0(e)X^(cr). 

We  compute  now  the  change  in   p{u  )   across   C;   by  the  mean 
value  theorem 

[p]  =  Pv;tw]  +  Pz^^^  • 
The  oscillation  of  w  and   z   is   0(e);   so  denoting  by   p 


w 
the  value  of  that  function  at   0,0  we  deduce  from  the  last 

three  relations  that 

(3.55)      [p^^^KcT)  =P^{x+(^^<7)  +x;(^^o')]  +0(£)x(^^cr) 

+  O(l)(x+^(or)  +  |X-(^^cr)|] 

By  definition,  an  approximate  characteristic  travels  either  with 
characteristic  speed  or  shock  speed.   In  either  case  we  have. 
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using  the  well  known  relation  that  shock  speed,  is  the  average 
of  the  characteristic  speeds  on  either  side,  see  [^  ], 

i^it)    =   p(Uj^(Xj,(t)+0,t))  -  I  p^Str;?'j^(t){l+0(e) 

We  recall  that  the  strength   (Str)   of  the  shock  on  "/C   was 
defined  as  the  absolute  value  of  the  change  in  w  across  it. 
We  can  now  compute,  using  (3-55)j  the  difference  in  the 
speeds  of  X-,,   and  'X^v   ^-'t  the  two  ends  of   cr,  : 

^2k(^^    -   ^Ik^^k^    =   Pw{^r~^^(^^    +x;(^^0')    -  |(Str  Xi+StrT^^) 
+   0(£)^^^^cr)    +   StrX]_    +  StrXgf 

+  0(l)(x;(^^c.)    +   (X-Z^^c)!]. 

We  let   k  tend  to   oo ;   all  quantities  tend  to  a  limit  for  all 
but  a  denumerable  set  of  t   and  we  get 

Lemma  J)  .6:      For  all  but  a  denumerable  set  of  t 

(3.56)   ■   X2(t)  -^(t*)  =  p^O(t)  +  0{l)fx+(o<t))  +  lx](o<t))|| 


where 


(5.57)     0(t)  =X+(a(t))  +x:(cr(t))  -^(StrXi(t)  +StrXo(t)) 


+  0(e)|x(att))  +  StrXi(t)  +  StrX2(t) 


We  turn  now  to  estimating  the  strength  of  shocks  on  charac- 
teristics.  We  say  that  a  characteristic  ^  which  is  the  limit 
of  a  sequence  of  characteristics  X.\^      runs  into  another  charac- 
teristic 0     which  is  the  limit  of  the  sequence  0,       for   "t  ^  t „ 
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if  all  but  a  finite  number  of  yC       run  into  0   before   time   t  ,. 
Let   I   be  an  interval  at  some  initial  time   t.,   containing 
0('(t.);   let  X-i  ,  Xp  t)e  a  pair  of  characteristics  issuing  from 
the  endpoints  of  1   which  do  not  run  into   X  for   t  ^  t„. 
Then  inequality  (2.30)  holds  for  infinitely  many  k;   passing  to 
the  limit  we  obtain 

(3.58)      StrX(tf)  1  lX~(l)|  +0(£)QCAj)- 

The  key  concepts  so  far,  such  as  characteristic  curve  and 
amount  of  wave  crossing  were  obtained  as  limits  of  the  approximate 
concepts  and  therefore  may  very  well  depend  on  the  particular 
subsequence  used.   It  may  happen  that  two  characteristic  curves 
"X   and   0   occupy  the  same  position  yet  they  do  not  run  into  each 
other;  there  may  even  be  a  finite  amount  of  wave  contained  in  the 
infinitely  thin  strip  between  themi 

We  show  now  that  nevertheless  the  quantities  in  question 
serve  as  upper  bounds.   According  to  the  discussion  at  the 
beginning  of  Section  2,  the  Riemann  invariant  w   increases  with 
increasing   x  across  a  right  rarefaction  wave,  decreases  or 
remains  constant  across  all  others.   Therefore  for  any  spacelike 
curve  o: 

Increasing  variation  of  w,   across   o-^^   =x   (c^v)» 
Since  increasing  variation  is  lower  semi-continuous,  we  conclude 
after  passing  to  the  limit  that 


79 


(3 '59)     Increasing  variation  of  w  along  cr  <_  IL    [cf)  . 

Let  X  ,  7(p,...  Xt^     "be  a  sequence  of  right  character- 
istics and  denote  by   cr.(t)   the  horizontal  Interval  between 

J 

"X.-,      and  y.   :   denote  by  L{t)   the  interval  between  "A-, 
and   Xm"   It  follows  from  the  definition  that  for  all  but  a 
countable  set  of  t 

X+(L(t))  =  YL   Xj:(0'j.(t)). 

We  shall  refer  to  this  property  as  the  additivity  of  X  . 

We  turn  now  to  the  conservation  laws,  i.e.  inequalities 
(3 '53)^  (3'5''+)«   These  were  derived  as  the  limits  of  approxi- 
mate conservation  laws  after  omitting  the  amount  of  wave 
cancelled  or  entering  the  boundary.   These  inequalities  can  be 
improved  by  replacing  Q,(/A)   by  Q,(A'  )   defined  as 

Q(A')  =  lim  Q^^^\). 

We  shall  refer  to  these  improved  inequalities  as  (3«53)'  and 
(3«5'^)';   the  important  property  of  Q(/\')   is  this:   Denote 
by  /\.(t)   the  strip  contained  between  '^.    ^      and  X-,-   and  under 
O'.(t),   by  TT('t)   the  strip  between  Xn   and  Xivr   and  under 

J 

o'{t).   Clearly  for  each  k 

so  letting  k  — >  oo   we  obtain 
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JZ  Q(/^i)  =  Q(Tr')  1  Q(TT) 


¥e  shall  call  this  the  addltivity  of 
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4 .   The  spreading  and  cancellation  of  rarefaction 
waves . 

In  Section  1  we  proved  that  if  X-,   and  Xn      ^^^ 
two,  say,  right  characteristics  of  a  dif f erent iahle 
solution,  then  X   and  Xp   diverge  at  the  rate   cW, 
where  W   is  the  increasing  variation  of  the  Riemann 
invariant   w  between  them  and   c   is  a  constant.  This 
estimate  was  used  to  prove  that  no  diff erentiable  space 
periodic  solutions  can  exist  for  all  time. 

In  this  section  we  derive  a  similar  estimate  for 
discontinuous  solutions,  stated  as  Theorem  4.1.   The 
quantity  W   is  the  amount  of  rarefaction  wave  between 
the  two  characteristics.   However,  unlike  the  case  of 
smooth  solutions,   ¥   can  decrease,  due  to  cancellation 
of  rarefaction  waves  and  shocks,  so  that  in  this  case 
the  estimate  on  the  rate  of  divergence  of  the  character- 
istics leads  not  to  a  nonexistence  theorem  but  to  an 
estimate  on  the  rate  at  which  W   decays  in  time. 

Just  as  in  Section  1,  we  shall  measure  the  distance 
between  right  characteristics  at  times  displaced  in 
the  direction  in  which  left  waves  propagate.   This  has 
the  effect  that  the  distance  between  the  two  right 
characteristics  remains  unaffected  (approximately)  as 
they  pass  through  left  waves. 

Let  /\  be  a  strip  bounded  on  the  left  and  right  by 
a  pair  of  right  characteristics  X   and  Xo*   Through 
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the  point  Xp(t),   t   we  construct  a  spacelike  line  segment 
a^      meeting  Xn   at  a  time   t"^ .   We  give   Oi   the  speed 
A-,  =  A{0,0)  -  0(e),  0(e)   large  enough  so  that  a,       is  space- 
like; the  choice  of   0(e)   determines   t* .   Let 


and 


Since 


and 


D(t)  =  X^i^)    -   ^i(t) 


D*(t)  =  yi^it)  -  y{^[t*) 


D*(t)  =  -  AQ(t*-t)(l+0(e) 


D(t)  ^   (p^-AQ)(t*-t)(l+0(e)), 


the  ratio 

(4.1)      D*(t)/D(t)  =  -^—  (l+0(e)) 

^o-Po 

is  approximately  constant.   Since  Xn   ^^'^  T^o   differ  by   0(e) 

D(t*)  =  D(t)  +  0(e)(t*-t) 
=  D(t)  +  0(e)D(t); 


this  implies  that 

(4.2)  D(t*)/D(t)  =  1  +  0(e). 

Following  the  proof  of  (1.12)  in  Section  1  we  obtain 

(4.3)  ^  D^(t)  =  {\{t*)-\[t))   —^—   (l+0(e)). 

^o-Po 

According  to  Lemma  3*6,  equation  (3-56) 
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(4.4) 


j(^{t*)    -   k^{t)    -  P^O(t)  +  0(l)(X+(a^)  +  |X^(o-^)  I). 


where  0(t)   is  defined  by  {3'57)«   In  order  to  estimate  the 
last  term  in  (4.4)  we  define 


t^  =  0(e)D*(t), 


with  the   0(e)   large  enough  so  that  the  line  segment   t, 
from  Xj^it*),    t*   to  X^{t    -  ^  t^) ,    t   -  ^   t^   is  timelike 


time  =  t* 


time  =  t 


Fig.  4.1 

Let   P"*,   be  the  region  bounded  by  cr.,    t,   and  ;(p.   Left 
waves  enter   [""*,  through   Xo  only  and  leave  through  c.      and 
T , .   Thus  using  (3 '54)' 


(4.5) 


+  , 


■,+/ 


|Xj(<^t^l  -  '^(^7)1  1  |Ej(pi  +  0(e)Q(rt) 

t+o 

|dX±(X2)|  +  0(e)Q(r^) 
'o 


t-t 
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We  define 

t       t 
(^.6)       at)  =  J  dX+  +  J  |dX~|  +  0(e)Q(/\^) 


o 


where  A  =  U  f"*  .   ¥e  can  write  (4.5)  as 
^   s<t   ^ 

t+0 


|Xj(<^t^'  -/   ^'^^^ 


t-t 
o 

We  substitute  the  above  estimate  and  (4.4)  in  (4.3) 
and  obtain 

(4.7)  ^  D*(t)  =  p^O(t)  ^  +0(1)  r    dC(t). 

^o-Po         t-t 

o 

This  equation  makes  sense  only  when  t-t  ^0,   which  is 

when  t^O(e)D*(t).   Since 

(4.8)  D*(t)  =  D*(0)  +  0(e)t, 

we  have   t-t   >  0   except  for 
o  —        -^ 

t  ^  t^  =  0(e)D*(0) . 

Using  (4.8)  again, 

D*(t^)/D*(0)  =  1  +  0(e). 

Since  D*   is  Lipschitz  continuous,  it  is  the  integral  of  its 
derivative  and  so  integrating  (4.7)  from   t^   to  T  we  obtain 

T  T   t+0 

(4.9)  D*(T)  =  (l+0(e))D*(0)  +  cj  0{t)    +  0{l)f   J       dC(s)dt, 

t^  t^  ^"^o 
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where 

A 
c  =  p  — ^—   (l+0(e). 
^  A  -p 

¥e  wish  to  eliminate   t^   from  the  single  integral  in  (4.9) 
To  accomplish  this  we  note: 

\J      Qdt|  ^  0(e)D^(0)sup^lO(t)  h 
o 

According  to  the  conservation  laws  {3.53)j  (3-54)   |o|   is 

bounded  by  the  total  variation  V^   of  the  Cauchy  data  over 

an  interval  of  length   O(t^);   so 


\J     Odt|  ^  0(£V^)D*(0) 


o 


Substituting  this  in  (4.9)  gives  us  the  equation 

(4.10)      D*(T)  =  (l+0(e)(l+V^))D*(0)  +cfo(t)dt 

T  t+o        o 

+  0{l)j  J  dC(s)dt. 

t„  t-t 
*  o 

The  double  integral  can  be  estimated  as  in  Section  1  by  a 

change  in  the  order  of  integration.   With  the   t-,(s)   introduced 

in  Section  1,  we  have 

T  J:+o  T_s+t^(s) 


J  J  dC(s)dt  ^j  j       dtdC(s). 


t^  t-t  o   s 


Using  the  estimate   t  (s)  =  0(e)D*(s)   from   Section  1,  we  bound 
the  integral  above  by 


0(e)J  D*(s)dC(s), 
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and  (4.10)  becomes 

T  T 

G{T)  -  r  D*(t)dm  <_  D*(T)  <  F(T)  +  rD*(t)dm, 
o  o 

where   m(t)  =0(e)C("t)   and 

T 
F(T)  =  (l+0(e)(l+V^)  )D^(0)  +  crO(t)dt, 

G(T)  =  (l-0(e)(l+V^))D*(0)  +  cj  0(t)dt. 

o 
Using  (1.28)  we  have 

(4.11)  G(T)  +  (l-e"^^^hF(T)  ^D*(T). 

Substituting  the  definition  of  F,   G   and   m  into  (4.11) 
gives  a  lower  bound  for  D*   in  terms  of  V^ ,  Cj  Ej  D*(0)   and 
O.   If  we  assume  that 

(4.12)  eV^  +  eC(T)   is  sufficiently  small 

we  can  simplify  that  inequality  to 

T 

(4.13)  c[l-small  no.]ro(t)dt  <_   D*(T). 

o 
Next  we  proceed  to  find  a  lower  bound  for  0{t);    using  the 

abbreviations 

xj(olt))  =  R(t),  X^(cr(t))  =  Sh(t) 

the  definition  (3-57)  of  0(t)   implies 

(4.14)  0(t)  >  [l-0(e)]R(t)  +  [l+0(e)]Sh(t) 

-  I  [l+0(e)][Str;^(t)+Str  Xglt)] 
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where   Str /(^  denotes  the  strength  of  the  shock  on  V. 

We  apply  now  the  law  of  conservation  of  right  shocks 
(3-53)  to  the  set  A(t)   bounded  by   c,  or  ,  X-,   and  Xp  • 

Dropping  the  amount  of  wave  cancelled  we  get 

(4.15)"    Sh(t)  >  Sh(0)  -  0(e)Q(/|'). 


+ 


Similarly,  the  conservation  law  (3 '53)   for  rarefaction  waves 
applied  to  the  set  bounded  by   cr,  ,  <y^,     X-i   and  Xo   yields 

(4.15)^    R(t)  >  R(T)  -  0(e)Q(/\'). 

Here  A    denotes  A(T). 

Finally  denote  by  Str  ^    the  maximum  strength  of  X 

max  X. 

and  A 2 •   Substituting  (4.15)—  into  (4.l4)  and  the  resulting 
inequality  into  (4.13)  yields,  after  division  by  T, 

(4.l6)     c[l-small  no.]R(T)  <_  ^*\^^    +  |Sh(0)I 

+  0(e)Q(A)  +  Str^^^. 

We  claim  now  that  inequality  (4.l6)  holds  without  the  term 

Str     on  the  right, 
max  ^ 

Proof;  We  draw  a  finite  number  of  right  characteristics 
X-,  =  ?-i  J  <...<  Iivr  =  "Kn      all  originating  in  the  interval 
[X-,(0)j  p(*^)]'   so  that  the  amount  of  rarefaction  wave  entering 
each  strip  /[  .      bounded  by   |  .,  ^  •  , -,   is  less  than   eQ,(/^' )  • 
Then,  according  to  the  law  of  conservation  of  rarefaction  waves 
the  amount  leaving  /I.   is  bounded  by   0(e)Q(/^'). 

J 


88 


Some  of  these  characteristics   ^  .   may  run  into   A-,   or 
Xp   still  in  ./{;      suppose  that   Ij.-i   does  but  ^^      doesn't 
run  into  %-,  ,   and  that   ^t,-,-]   does  but   ^x.   doesn't  run  into 
Xp •   According  to  the  additivity  of  the  amount  of  rarefaction 
wave  which  crosses   o'(T)   between   ^f-,   and   7(^  ^  R(T)   is  the 
sum  of  the  amount  crossing  between   X-i  ^-^^      ^j>      ij     ^■'^'^ 
^    ,      and  ^y   and  Xn-      Since  Xn  and  |j_-,  have  coalesced  on  (y{T) , 
the  first  amount  is  equal  to  what  is  leaving  A-      and  this  is 
0(e)Q,(/\');   the  same  estimate  holds  for  the  third  amount. 

'J 


Thus   R   (T),   the  amount  crossing  cf{l)      between   | -.   and 


i        differs  by   0(e)Q(/\')   from  R(T): 

(4.17)  R(T)  <  R'^^{T)  +0(e)Q(A'). 

¥e  apply  now  (3.58)  to   |j   and  to  i    ,      and  the  strip 
-A;   we  obtain  that  the  strength  of  the  shock  on  each  is  bounded 
by   |Sh(0)  I  +  0(e)Q(/l): 

(4.18)  Str^^^  <  lSh(0)l  +  0(e)Q(yi'). 

¥e  apply  now  inequality  (4.l6)  to  the  strip  bounded  by   |^   and 


J 


I   and  obtain,  using  (4.l8),  that 


c[l-small  no.]R'^^(T)  ^  D-^  (T)/T  +  |Sh(0)|  +0(e)Q(/i'). 
Using  now  (4.17)  and  the  fact  that  B*-^     <_  D^   we  obtain 

(4.19)     c[l-small  no.]R(T)  ^  5!_m  +  sh(0)  +0(e)Q(y\'). 

T 

In  the  above  derivation  we  have  assumed  that  (4.12)  holds 
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If  all  or  all  but  one  of  the   £.   run  into   |-^   or  ^^, 
then  R(T)  ^0(£)Q(/\')   follows  from  {3'53)'^- 

We  are  now  ready  to  prove  the  main  result  of  this  section: 
Theorem  4.1:   Let  X^      and  Xp   be  a.  pair  of  right  charac- 
teristics,  TT(T)   the  strip  bounded  by  Xi      ^'^'^    Xp^  0  ^  t  <_  T. 
Denote  by  L ( t )   the  interval  [X-.  ( t )  ,Xn  ( 't )  ]  ^   its  length  by 
|L(t)|.   Denote  by  V  (TT)   the  total  variation  of  the  Cauchy 
data  over  the  domain  of  determinacy  of  7T(T).   Suppose  that 

(4.20)  eV  (TT)   is  sufficently  small; 

then 

(4.21)  ax;;;(L(T))  <  Mm  +  o(£)Q(TT)- 

^  T 

The  constant   a  can  be  taken  arbitrarily  close  to   [l+0(e)]p  (O) 
by  making  the  bound  in  (4.20)  small  enough. 

Note:   Clearly  this  inequality  asserts  that  either 
characteristics  curves  spread  apart,  or  the  amount  of  rarefaction 
wave  between  them  decays. 

Proof:   ¥e  divide  TT  into  alternate  strips  of  predominantly 
shock,  respectively  rarefaction  flow.   More  precisely,  we  divide 
the  initial  interval  L(0)   into  a  finite  number  of  subintervals 
L.(0)   so  that   |L.(0)|  <  e   and  L.(0)   contains  less  than 

J  J  J 

6  .   amount  of  right  shock,  respectively  rarefaction  wave,  depend- 

J 

ing  on  whether   j   is  even  or  odd.   If  the  initial  data  are  of 
bounded  variation,  we  can  make  the  subdivision  so  fine  that 

>~  5  .  =  6 
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is  less  than  any  preassigned  number. 

Denote  "by   2,  .   the  characteristic  issuing  from  the  j 
subdivision  point,  and  by   /I.   the  strip  bounded  by   ^., 
^    ,  L.(0),   and  cr.{T.),    T.   so  chosen  that   T^  =  T.   For   j 
odd, 

x;(Lj(0))  <  5j; 

then  according  to  the  law  of  conservation  of  rarefaction 
waves  (3  -53  )  ' 

(4.22)  .  X^(L.(T))  ^  5  .  +  0(e)Q{/\;-)- 

For   j   even  we  apply  inequality  (4.l9);  using  the  fact  that 

|Sh(0)l  <  5.   and  that  according  to  the  limiting  form  of 

Theorem  2.1,   C(T)   is  bounded  by  V  (/\)   we  get  from  (4.20) 

that  (4.12)  is  satisfied,  so 

D^(T.) 

(4.23)  •    c(l-small  no.  )R  .(T  .)  ^  -^ ^  +  5  .  +  0(e)Q'. . 

J  J    J  m  J  J 

¥e  apply  now  the  conservation  law  for  rarefaction  waves  to  the 
region   p*.   shown  in  figure  4.2: 

J 


Fig.  4.2 
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We  eet 


xJ(Lj(T))  <  R.{T.)    +  0(e)Q(  P^. ) 


According  to  (4.1) 

A 
D*.(T  )  -  iL  (T)l  ^-^  (l+0(e)). 
J   J      J     ^o  ^-^o 

It  is  easy  to  see  that  since   |L.(0)|  <  e, 

J 

T  .  =  T{l+0(e)  ). 
Substituting  these  three  relations  into  (4.23)  •  and  recalling 

J 

that 

A 
_o 

■o~''o 

we  get 


C  =  p,^  —^^—     (l+0(£)) 


|L.(T)|  _   I 

(4.24)      p  (1-small  no.)X+(L.(T)  )  <  — ^ +  B  •  +  0(e)[Q(AJ^^(  FT)] 

J"*  ■'•J  rpj  Jd 


We  add  up  the  inequalities  (4.22)  .   and   (4.24)  .;  by  additivity  of 

J  J 

amount  of  wave  crossing 

xJ:(Lg(T))^  =  X+(L(T)) 

by  additivity  of  length 

i:  lL.(T)l  ^  lL(T)|; 
J 

by  construction 

5  .  =  B 
J 
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Finally,  according  to  the  additivity  of 


Q(/\^-)  +  ki   r^.)  1  Q(TD 


Substituting  these  relations  into  ^  (4.24)  .   we  obtain  that 

J 

p  (l-small  no.)X^(L(T))  <  ^!^   +  5  +  0(e)Q(TT). 
'  w  r        —   i 

Since   6   is  arbitrary,  inequality  (4.21)  of  Theorem  4.1  follows. 

¥e  claim  that  inequality  (4.21)  holds  for  any  interval 
L(T),   not  necessarily  bounded  by  two  characteristics;  for 
this  purpose  we  construct  two  characteristics  X,   and  Xn      ^^ 
that   L(T)   lies  between  /(-,{T)      and  Xp(T).   Since  the  speed  of 
propagation  of  any  right  characteristic  is   p(0)  +  0(e),   we  can 
choose  /t-,   and  Xn      so  that   ^.(T)   and  Xp(T)   lie  within  a 
distance  of  0(e)T   of  L(T).   Denote  by  TT(L)   the  strip 
contained  between  X-,   and  Tfo   so  chosen.   Next  we  draw  a  finite 
number  of  characteristics  X-,  =  I-,  <  ^^  ^•••'^  ^-.^   =    '^ n      ^^'^ 
originating  in  the  interval   [XnlO),  Xp(0)]   so  that  the 
amount  of  rarefaction  wave  entering  each  strip  bounded  by   ^ ., 
4  .     is  less  than   eQ,(TT)  •   Denote  by   J  and   K  the  indices 

J  "T  X 

such  that   ?j_-|(T)   does  not  lie  in  L(T)   but   |.(T)   does, 
and  similarly  so  that   ^|^(T)   lies  in  L(T)   but   ?k+i(T)  doesn't. 
Thus,  since   L(T)   is  contained  in   ^j.i^^)^  ^K+1^^^'   ^®  have  by 
additivity  of  waves  crossing  that  at   t  =  T 

(4.25)    x>)ixj(4j_^,|^^^)  =x+(|j_^,ej)  +x^(4j,|^)<(|^,Ik^-,), 

Applying  the  conservation  law  (3-53)  we  get  using  the  bound  imposed 
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on  X'^ii  .,i  .  ,.)      at   t  =  0   that  at   t  =  T 


(4.26)      ^r^^J-l'^J^  +^r(^K'^K+l^  lO(e)Q(rD 


We  apply  now  Theorem  4.1  to  the  strip  "bounded  by   |^,  |„; 
we  get  an  estimate  for  X  (|-|-,^„)   in  terms  of   |  ^  -^  |   which 
is   <_  L.   Combining  this  inequality  with  (4.25)  and  (4.26)  we 
obtain  inequality  (4.21)  for  an  arbitrary  interval.   In  light 
of  the  observation  (3-59)  we  can  express  this  result  as  follows; 

Theorem  4.2;   Let   L  be  any  x-interval  at  time  T;   let 
"pr  be  any  characteristic  strip  containing  L,   and   suppose 
that 

(4.27)  eV  (TT)   is  sufficiently  small. 

Then 

(4.28)  p  (1-small  no.)   Increasing  variation  of  w   over  L 


w 


<  M  +  0(e)Q(TT) 


For  purposes  of  easy  reference  put  (4.28)  into  the 
following  form: 

IlI  -r-r 

(4.29)  Increasing  variation  of  w   over  L  ^  K  -• — L  +  K-|eQ(n). 

O  m         -'- 

The  total  variation  of  w   is  at  most  twice  its  increasing 
variation  plus  its  oscillation;  since  the  total  oscillation 
of  U   is  E,      that  of  w   is   0(e).   We  get  a  similar  estimate 
for  the  total  variation  of   z.   Combining  the  total  variation 
of  U  =  fu,v}  : 

(4.30)  Total  Variation  of  U^^^  over   L  <_  K  -L^  +  k^e  +  Ko£Q(TT) 

2  rp        2        2 
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5 .   The  existence  and  decay  of  solutions  with  arbitrary 
Cauchy  data  near  a  constant  state. 

In  this  section  we  solve  the  initial  value  problem  with 
initial  data  that  are  merely  measurable,  and  which  differ 
little  from  a  constant  in  the   L   norm.   We  show  that  even 
if  the  initial  data  have  infinite  total  variation,  the 
corresponding  solution  U(x,t)  ^  j  u{x,t)  ,v{x,t)j    has  for  every 
positive   t   finite  variation  which  decays  in  time. 

For  single  equations  this  decay  of  total  variation  has 
been  proved  by  Oleinik,  [  6  ] ;  the  decay  is  a  quantitative 
measure  of  the  smoothing  which  the  initial  data  undergo,  see 

Throughout  this  section  C,  K  and  H,   with  various 
subscripts,  denote  constants  which  depend  only  on  the  equation 
under  consideration,  and  not  the  solution  in  question. 

Theorem  5.1:   We  are  given  bounded,  measurable  initial 
data  U(0);   denote  its  maximum  by  r\    : 

(5.1)  I|U(0)|1^  =  7^2. 

If  T]   is  small  enough,  there  exists  a  solution  U  =  ^u,v? 
of  the  system  of  conservation  equations  (1.1  )  with  initial 
values  U(0);   this  solution  satisfies  the  decay  law 

H  |l| 

(5.2)  Total  Variation  of  U(x,t)  over  L  <  — = 

~   t 

for  every  x-interval  L  whose  length   11]   is_  ^  t .   Furthermore 


(50)  |u(x,t)  1  <  H2TI 


95 


We  shall  rely  on  the  following  existence  theorem  and 
estimates  derived  in  [ 2  ]  : 

Existence  Theorem;   Suppose  that  initial  data  U(0) 
are  prescrioed  subject  to  the  following  restrictions; 

i)   1|U(0)||^  1  C^ 
(5.4) 

ii)   ||u(0)|l^rTotal  Variation  of  U(0)j^  ^  C^. 

Then  there  exists  a  solution  U(x,t)   with  initial  value  U{0) 
which  satisfies 

i)   Total  Variation  of  U(t)  1  C^  Total  Variation  of  U(0) 

(5.5) 

ii)   |u(x,t)|  <  C^||U(0)||^, 

provided  that   C-,  ,  Cp   are  sufficiently  small  and   C^,  C^, 
sufficiently  large. 

It  is  convenient  to  assxime  for  what  follows  and  through- 
out this  section  that  all  characteristic  speeds  are   <  1. 
This  permits  us  to  set   r,  =  s,  in  the  construction  of  approxi- 
mate solutions. 

Local  Existence  Theorem:   Suppose  the  initial  data  U(0) 
satisfy  condition  (5.^)-^  b'U-t  (5.'^)--  only  over  sub-intervals 
of  the  X-axis  of  length  5.   Then  a  solution  U(x,t)   with 
initial  values   U(0)   exists  for   t  <  6/2;   (5-5)  •   is  satisfied 
over  intervals  of  length   |l|  <  6  -  2t. 


9^ 


Sketch  of  proof:   Let   J  be  a  subinterval  of  the  x 
axis  of  length   6;   define   U-j.(0)   to  be  equal   U(0)   on   J, 
constant  outside  of  J.   These  initia.l  data  satisfy  (5-'^); 
therefore  a  solution  U-|.(x,t)   exists.   Suppose  that  in  the 
approximation  procedure  used  to  construct  these  solutions  we 
employ  the  same  sequence  of  time   and  space  increments  and 
the  same  random  numbers   a  .   Since  the  propagation  of  influ- 
ence  in  the  approximate  solutions  is  —  =  1,   it  follows  that 
if   J-,   and   Jp   are  overlapping  intervals,  the  corresponding 
approximate  solutions  are  equal  in  the  domain  of  determinacy 
of  the  overlap.   The  same  is  true  for  the  limits  if  we  choose 
the  same  subsequence;  this  is  possible  if  we  employ  only  a 
denumerable  set  of  intervals   J.   We  set  now 


U(x,t)  =  Uj(x,t)   in  the  domain  of  determinacy  of 
J.   This  completes  the  proof  of  the  local  existence  theorem. 

We  break  up  the  proof  of  Theorem  ^.1   into  four  steps. 

Step  1:   It  suffices  to  prove  the  theorem  for  solutions 
whose  initial  data  are  uniformly  smooth. 

Reason:   Arbitrary  bounded  initial  data.  U(0)   can  be 
approximated  by  a  sequence  of  mollified  data  U  .   defined  by 

U^(0)  =  U(0)^j^; 

here   j  (x)  =  nj(nx),  j   some  smooth  nonnegative  function  with 


This  argument  is  necessary  since  we  sadly  lack  a  local 
uniqueness  theorem. 
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compact  support  and  total  integral  1.   Clearly 


i)   llUn(0)||^  1  l|U(0: 


il)   U^(x,0)  -^  U(x,0)   a.e 


ill)   There  is  a  number  M^   such  that  \-^-   U^(x,0)l  <  M  . 

Suppose  we  have  shown  that  the  solutions   U    satisfy 
inequalities  (5.2)  and  (SO);  then  by  Kelly's  theorem  we  can 
select  a  subsequence  which  converges  to  a  limit  U;   this  limit 
has  initial  values   U(x,0)   and  satisfies  inequalities  (5'2) 
and  (5.3). 

Step  2:   It  is  enough  to  prove  (5-2)  for  intervals   L 
of  length  equal  to   t,  with  H-, /2   in  place  of  H-,  . 

Reason:   If   (5.2)  holds  for  intervals  of  length   t,   it 
also  holds  for  intervals  of  length  mt,   m  an  integer.   It 
follows  then  that  (5.2)  holds  for  intervals  of  length  mt, 
m  >  1  with  H-,   in  place  of  H-, /2 . 

Step  3 •   Local  existence  of  the  solution. 

Suppose  we  are  given  uniformly  smooth  data  U(0);   we 
choose   6   so  that  the  total  variation  of  U(0)   over  intervals 
of  length   &   is   <_  1.   We  restrict   t]   to  be  so  small  that 

Then  the  hypotheses  (5.4)  of  the  local  existence  theorem  are 
satisfied  over  interval  of  length  35.   The  local  solution  exists 
in  the  strip   t  _<  ^  ,   and  satisfies  inequalities  (5*5)  in  each 
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right  triangle  whose  hypotenuse  lies  on  the  x  axis  and  is 
35   long.   Clearly  property  (5*5)  ••  implies  property  (So) 
provided  that   -q   is  so  small  that 

{5.5).  asserts  that  the  total  variation  of  U(t)  along  an  interval 
of  length   (35-2t)   is  ^  ^^>      this  implies  (5.2)  with   |l|  =  t 
for  t  ^  5,   provided  that 

(5-6)  C^  <  H^/2. 

This  represents  a  limitation  on  H-,  . 

Step  4:   Inductive  proof  that  a  solution  U(x,t)   exists 
for   t  ^  2  5   and  satisfies  (5  •2)  and  (5-3)  there. 

Step  3  shows  that  the  assertion  is  true  for  n  =  1. 
Suppose  now  that  (5-2)  and  (50)  hold  for   t  ^  2  5.   Then  the 
data  U(x,2  5)   satisfy  conditions  (5-'^)  over  intervals  of 
length   2    &,   provided  that   r|   is  so  small  that 

H^Ti  1  C^    ,      H2il(2H^)^  <  C^. 

If  so,  according  to  the  local  existence  theorem  the  solution 
U(x,t)   exists  in  the  strip   2  6  ^  t  ^  2    5.   What  we  have 
to  show  is  that  (5*2)  and  (5*3)  are  satisfied  there. 

We  note  first  that  it  follows  from  the  local  existence 

theorem  that  (5*5)  ••  holds;   combining  this  with  (5-3)  at 

n 
t  =  2  5   we  conclude  that 

(5-7)       IIuHqo  1  C^Hg-n   for   t  <  2'^+^&. 
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¥e  turn  now  to  estimating  the  total  variation  of  U 
for  t  _<  2    5.   According  to  inequality  (4O0),  applied  to 
an  interval  L   of  length  t,    with  initial  time   t/2   and 
T  =  t/2, 


(5.8) 


Total  Variation  U(t)   over  L  < 


<  K. 


L 


t/2 


+  K^e  +  K^eC 


=  2K2  +  K2(14^)e, 


where  Q   is  the   Q- measure  of  "ff  ,      the  domain  of  dependence 
of  L   in  the  strip   (t/2,t): 


L 


t/2  .- 


> 


Fig.  5.1 


According  to  the  estimate  (2.l4) 


(5.9) 


p 


V  the  total  variation  of  U   over  the  base  of  "|T-   Since  the 
base  of  "fT  has  length  2t   and  is  located  at  time   t/2  <  2^6, 
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according  to  induction  hypothesis  (5.2)  is  applicable  and  gives 


(5-10)  V  ^  4h^. 

Combining  this  with^5.9)we  get 

(5-11)  Q  1  32H-^^; 

substituting  this  into  (5«8)  as  well  as  the  estimate  (5-7) 
for   e  we  get 

(5.12)     Total  Variation  U(t)   over  L  <_  2K^   + 
+  K2(l+32H^^)C^H2Ti. 

For  T]   small  enough  the  second  term  on  the  right  in  (5-12) 
is  _<  Kpj   for  such   t] 

Total  Variation  U(t)   over  L  <_   3Kp. 

Clearly  with  H-,  =  6Kp   the  right  side  is   1/2H-,   which  accord- 
ing to  step  2  implies  (5-2)  for  all  L   of  length  >_   t. 

We  note  that  for  the  validity  of  estimate  (4.30)  it  was 
necessary  that  vl|u||    be  small  enough.   Using  estimates  (5-7) 
and  (5 '10)  we  see  that  the  above  quantity  is  <_  hY{^Q,uR^T\,    which 
is  small  enough  if  t\      is.   This  completes  the  proof  of  (5 '2) 
for   t  <_  2    6 . 

We  turn  now  to  proving  (5-3)  for  t  <_  2    5;   this  calls 
for  estimating  ||u(t)l|   ,  or,  what  is  the  same,  estimating 
||w(t)||    and   l|z(t)ll   .   Inequality  (4.29)  gives  an  estimate 
for  the  increasing  variation  of  w: 
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Increasing  variation  of  w  on  L 


L 


(5.13)    1  Kq  ^^  +  k-l£Q(TT). 

where  JJ  is  the  domain  of  dependence  of  L   in   (t-T,t). 
(5-13)  alone  is  not  sufficient  to  estimate   11  w  11   ,   unless 
combined  with  some  absolute  estimate  for  w.   Such  an  absolute 
estimate  is  provided  by 

Lemma  5 » 1 •   Let   L  be  any  interval  located  at  time 
t  :   then 

(5.1^)     /  |w(x,t)  |dx  <  K2,(  lLl+t^)|iu(0)II^, 

'l 

Kn   independent  of  U;   a  similar  estimate  holds  for   z. 

We  postpone  the  proof  of  Lemma  5'1  and  show  first  how 

to  use  it  to  estimate  w.   Let  x  ,  t   be  any  point  and  denote 

w(x  ,t  )   by  w  :   suppose,  say,  that  w    is  positive.   Let 

L  be  an  interval  at  time   t   whose  right  endpoint  is   x  ,  t  . 

o  to      jr-  00 

Using  (5 -15)  we  get  the  following  lower  bound  for  w   on  L: 

^0-^0       K^£Q(TT)  1  w 

Integrating  this  inequality  over  L  and  using  (5'1^)  we  get, 
after  dividing  by   |l|,   that 

We  recall  now  from  (5.I)   that   i|u(0)||    =  t)  ;   choosing 
|l|  =  Tit   we  get  from  the  above  inequality  that 
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(5.15) 


w. 


<    K^Ti       +   K^^T)    +   K^   ^  Ti    +   eK-j_Q(TT) 


¥e   recall  that      e      is   the  maxlmiiin  of 


U 


in  IT  and  there- 


fore is  certainly  >_  constant  multiple  of  w  .   Therefore 
(5 '15)  cannot  furnish  an  upper  bound  for  w   unless  the  coef- 
ficient of   e   on  the  right  in  (5-15)  is  small  enough,  i.e. 
unless  Q,(TT)   is  small  enough.   Our  only  available  estimate 
for  Q  is   (5 '9);  and  our  only  available  estimate  for  V   is 
(5-2);  therefore  our  best  estimate  for  Q   is   2H-,    which 
certainly  isn't  small.   We  have  to  search  for  special  domains 
with  small  Q-measure;  we  proceed  as  follows: 

Let   y,  s   be  any  point  in  the  strip   2^6  <  s  <  2^  &; 
denote  by   p*  the  rectangle: 


y-s<x<y  +  s,        2^''"6<t<   2^5 


2^+15  <f 


2^5 


y-s 


y+s 


Fig.  5.2 
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The  length  of  \~'     in  the  x   direction  is   <_  2.2    5;   so 
according  to  (5 '2)  the  total  variation  V   of  U  along  the 
lower  edge  of   [""  is   <_  8H-,   and  hence  by  (5-9) 

(5.16)  Q(  p)  1   128H^^. 

We  divide   p""  hy  horizontal  lines  into  N   congruent  rectangles; 
each  subrectangle  has  vertical  width  2  ~  5/N;   at  least  one  of 
them,  call  it   p.,,   must  have   Q-measure  <  Q(  r~')/N;   by  (5.16) 


N 


2 


5-17)        Q(  r^)  1  i28h^Vn. 


Denote  by  t   the  time  on  the  upper  edge  of   pL   and 

let  X  ,  t   be  any  point  on  that  edge  which  lies  in  the  domain 

of  dependence  of  y,  s .   Let   L  be  an  interval  at  time   t 

whose  right  or  left  endpoint  is   x   and  whose  length  is   r^t  ; 

clearly,  for  r\   <    l/4   the  domain  of  dependence  XT  of  L   in 

ervax  y.    -x,.  ;,   x  =^^-^^ 
^    o      '    o   ' 

So  by  (5.17) 


the  time  interval   (t^-T,t^),   T  =  2^^  b/N,   belongs  to   p  . 


Q(TT)  1  Q(  r^)  1  128h^2/n. 

Denote  as  before  w(x  ,t  )   by  w  ;   according  to  (5 "15) 
with  T  =  2    6/n  and  the  above  estimate  for  QlTT) 

(5.18)      Wq  ^  K^^Ti^  +  K^ri  +  2K^Nti  +  eK-^  128h^^/N. 

A  similar  estimate  holds  for   z   =  z(x  ,t  )   at  any  point  x  ,  t 

o      o'  o        "^  ^       o'   o 

in  the  domain  of  dependence  of  y,  s.   Since   u  and   v   are 

We  have  to  avoid  having  points  with  positive  Q-mass  on  the 
horizontal  dividing  lines. 
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Llpschitz  continuous  functions  of  w  .and   z   we  get  a  similar 

estimate  for  U   =  U(x  ,t  ): 

2 
eKH, 

|U  I  _!  KNt]  +  ^  , 

"  N 

where   K   is  some  constant  concocted  of  the   K.,  and  where  we 

J 

n 
have  used  that   t]  ^  1?  N  ^  1.   Since   t   <  2  B,   "by  induction 

hypotheses  (5-3)  holds  and  can  he  used  as  estimate  for   e.   So 

KH^HpTi 

lu    I    _^  KNri    + ^    . 

N 


Clearly  if  we  choose 


and  set 


N  =  2KH^^C2^ 


Hg  =  2KN 


then  the  above  inequality  gives 
(5.19)      |U J  <  -^-  . 

Denote  by  1^   the  domain  of  dependence  of  y,  s   at   t  =  t  . 
By  induction  hypothesis  (5 '2)  and  {5-3)  hold  there;  it  follows 
that  for  T]   small  enough  the  hypotheses  (5''^)  of  the  local 
existence  theorem  hold,  and  therefore  so  does  the  conclusion 
(5 '5)  in  the  domain  of  determinacy  of  I^.   Since   y,  s  belongs 
to  that  domain  of  determinacy  we  conclude  from  (5.5).-  and  (5'19) 


that 


|U(y,s)|  <  Cj|u(t^)||^  1  C^  -^  =  HgT^, 

C4 
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Since   y,  s   is  an  arbitrary  point  in  the  strip   2  6  <  s  <  2    5, 

this  completes  the  proof  of  (5«3)« 

We  turn  now  to  the  proof  of  Lemma  ^.1;    we  need 

Lemma  5 ♦ 2 :   For  any  solution  U   the  maximum  of  w  and 

the  minimum  of   z   is  reached  at   t  =  0. 

Proof:   Let  U.   and  U   denote  the  states  on  the  left 

and  right  respectively  in  a  Riemann  initial  value  problem,  and 

let   U   denote  the  intermediate  state,  which  is  connected  to 
m 

U.   across  a  left  wave,  and  to  which  U   is  connected  across 
i  r 

a  right  wave.   The  Riemann  invariant  w  doesn't  change  across 
a  left  rarefaction  wave  and,  according  to  ass\imption  (2.4)  in 
Section  2,  decreases  across  a  left  shock.   Thus 

w(Ujj^)  <  w(U_g). 

Similarly  z   doesn't  change  across  a  right  rarefaction  wave 
and  according  to  (2.4)  decreases  across  a  right  shock.   Thus 

z(U  )  >  z(U  )  . 
^  m  —   ^  r 

Clearly  the  same  inequalities  hold  for  any  state  U'   in  the 
rarefaction  waves,  if  any.   This  proves  Lemma  5*2  for  the 
Riemann  initial  value  problem.   Since  our  approximate  solutions 
were  pieced  together  from  solutions  of  Riemann  problems.  Lemma.  5*2 
holds  for  all  approximate  solutions,  and  so  for  their  limits  as  well. 

Remark:  If  the  sign  in  any  of  the  inequalities  (2.4)  is 
changed,  the  corresponding  phase  in  Lemma  5*2  has  to  be  changed 
from  maximum  to  minimum,  or  conversely. 
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In  what   follows   we   shall   employ   the   abbreviation 


(5.20)  llU(0)lloo    =  ^' 


Since  w  and   z   are  smooth  functions  of  U,   for  any  x 

w(x,0)  <  k|u(x,0) I  ^  KM 
and 

z(x,0)  >_  -   k1u(x,0)|  >  -  KM. 

Using  Lemma  ^.2   we  deduce  this 

Corollary:   For  any  x  and  any  t  >  0 

(5.21)     w  ^  KM,   z  1  -  KM. 

This  corollary  can  be  stated  geometrically:  the  points 
(w,z)   lie  in  the  second  quadrant  of  the  w,  z  plane  with 
respect  to  the  point   (KM, -KM).   Since  the  mapping 

(5-22)     w,z  ->  u,  V 

is  smooth,  it  maps  this  quadrant  into  some  angle  with  vertex 
at.  the  image  of   (KM, -KM).   Replacing  if  necessary  the  variables 
u,  V  by  suitable  linear  combinations  we  can  achieve  that  the 
angle  to  which  u,  v   is  confined  is  the  first  quadrant.   The 
analytical  form  of  this  statement  is: 

(5-23)      -  KM  ^  u,   -  KM  <  v; 

here  we  have  used  the  fact  that  (5.22)  maps  the  origin  into 
the  origin,  and  we  have  enlarged  the  constant   K. 
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It  follows  immediately  from  (5 -23)  that 


(5.24) 


u|  <  u  +  2KM,  |v|  <  V  +  2KM. 


We  take  now  the  original  conservation  equations  (l.l) 


u^  +  f^  =  0,   v^  +  g^  ^  0, 


add  them  and  integrate  them  over  the  trapezoid  shown  on 
figure  5 '3: 


/N 


o 


Fig.  5  0 


Integration  by  parts  gives 


(5-25)     f  {u+v)dx  =  j    (u+v)dx  -  j    [ ( u+v) -( f+g) t^] 


dx. 


o 


Since  we  may  add  an  arbitrary  constant  to  f   and  g,,      we  can 
assume  that 

f(0,0)  =  g(0,0)  =  0. 
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Since   f  and   g  depend  smoothly  on  .u   and  v, 

|f  I  -I  lg|  ^  K(  |u|  +  |v|  ). 
Using  this  estimate  and  (5.24)  we  get 
(5.26)      (u+v)  -  (f+g)t^  >  ( lu|+lv| )(1-K|t^| )  -  2KM 

We  choose  now  the  sides  of  the  trapezoid  so  that   1  =  Kit  |; 

this  makes  the  length  of  each  side  in  the   x   direction  =  Kt  . 

^  o 

Using  (5.26)  in  the  second  integral  on  the  right  in  (5-25)? 
and  (5 '24)  on  the  left  we  obtain 


j{   I  u  I  +  I  V  I  -4km  )  dx  <_  M  1  u  I  +  I V  I  )  dx  +  4KMKt^ . 
L  L 


1  .6 


/(|u|+|v|)dx  <(4kL  +  2Kt^  +  4k  t^lM. 
L 

Since  w  and   z   are  smooth  functions  of  u,  v,   and  since 

M   is  an  abbreviation  for   ||U(0)|1   ,  inequality  ( 5 .  l4 )  of 

Lemma  5-1  follows  from  the  above  inequality.   This  completes 

the  proof  of  Theorem  5-1- 

We  turn  now  to  the  study  of  periodic  solutions: 

Theorem  5-2:   Suppose  that  the  initial  values  of  U  are 

periodic  with  respect  to  x-,   with  period   p.   Then 

a )  There  exists  a  solution  U ( x , t )   which  is  periodic 
with  respect  to  x,   with  period  p,   for  all  t. 

b )  The  average  U   of  U   over  a  space  period  is  inde- 
pendent of   t  . 
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c)   V(t),   the  total  variation  of  U   over  a  space  period 
satisfies  the  decay  law 


(5.27) 


V(t)  ^  %  , 


-,    "iP 
d)   |u(x,t)  -  U|  <  -^—    . 

t 

Proof:   It  follows  from  the  construction  described  In 
Section  2  that  if  the  initial  value  has  space  period  p,   so 
do  all  approximate  solutions  for  all  t,      and  therefore  so  does 
their  limit.   This  proves  a). 

To  prove  b)  Integrate  the  conservation  law  for  U   over 
the  rectangle   (0,p)  x  (0,t)   and  observe  that  because  of 
periodicity  the  flux  integrals  cancel. 

To  prove  c)  we  take  the  Interval  L  in  Inequality 

(5.1)  to  have  length  mp,   m   so  large  that  mp  ^  t .   The 

total  variation  of  U   over  such  an  L   is   mV(t),   so  by   (5-1) 

H-,mp 
mV(t)  ^  — — 
t 

which  Implies  (5-27)- 

Part  d)  follows  from  b)  and  c)  since  total  oscillation 
is  bounded  by  the  total  variation  per  period. 
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method,  or  process  disclosed  in  this 
report. 
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